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PREFACE 


The  sample  design  and  passenger  flow  estimation  work 
described  in  this  study  was  performed  under  Project  Plan 
Agreement  UM-937,  sponsored  by  the  Urban  Mass  Transporta- 
tion Administration,  Office  of  Planning  Management  and 
Demonstrations,  Office  of  Transportation  Management.  The 
work  was  undertaken  to  develop  a methodology  of  estimating 
the  distribution  of  passenger  flow  for  the  Light  Rail 
Transit  at  the  underground  stations.  The  techniques  are 
developed  in  the  context  of  the  Massachusetts  Bay  Transit 
Authority  - Riverside  Line,  but  could  be  expanded  to 
accommodate  other  transit  properties  with  appropriate 
modifications.  TSC  Technical  Program  Manager  is  Mary 
Roos , DTS-722. 


i i i 


METRIC  CONVERSION  FACTORS 


i 1 

I ^ 

a 

: : * 

t I s 

■ ■ tu 

* a -• 


1 f c l i 


Ihli 


8 ♦ ^ 
o o M » o 


~ \ 5 

I - ‘ * I 

i $ { { I 


"i  Vi 


2 • 


z » » a~*  i 


I i .1! 

ill!  ill  Hlili 


1 


S ^ 8 S « 

® “ a - 


j 1 1 f 

nil  ili  liiiii 


I s 

f f 


( 

. 1 


Is 


]l 


I S i f S '^''f  1 3 


■ — - f i 


-s 


ct 

llllllHI 

tt 

lllllllll 

u 

Illllilll 

ot 

llllllltl 

(1 

Illllilll 

tl 

llttlllll 

41 

ijilll 

91 

IllJlIll 

11 

Illllilll 

Illllilll 

Cl 

Illllilll 

Cl 

llltlllll 

It 

lllllllll 

ot 

imliiii 

• 

iinLi 

i 

iiiiLi 

4 

1 

9 

lllllllll 

9 

lllllllll 

* 

lllllllll 

C 

iiiiLii 

c 

ijllll 

lllllllll 

•lyr 

T|T 

T|T 

T|T 

•|'|T 

T|*l' 

T|T 

T|T 

M'|T 

T|T 

T|T 

T|T 

T|T 

T|T 

T|T 

T| 

• 

7 

t 

9 

4 

3 

a 

t 

•Kll 

S S * S l-i-vl  1 


lEE EC 


I i 

f I 
i i 


! 

• 

1 • 

i f 

» 

it 

t i . 
_ _ • 

j. 

. f - j ? 

ill  7 ^ 

• 

1 

a 

3 \ 

11 

n 

111 

ill 

iinfllll 

i 1 

• 

(J 

s a 


^ a o - 2 

s 

•«  8 • « 9 

« oi  b ««  6 ^ 

4 

• • s 

R ® ® S 

w 

• 

M 

• 

^ r * • S * " 

K 

< 

m 

e 

< 

y% 

< 

a 

III! 

> 
, i 

m 

X 

a 

35  ^ 

Ij!I  llni  1 


flJHiiil 


if 

:!n 

s 


s « T I 


i"«Vi 


3 » 


ill » u a r"\ 


IV 


TABLE  OF  CONTENTS 


Section  Page 

I.  INTRODUCTION 1 

II.  THE  SURVEY  PLAN 2 

III.  DATA  EDITING 8 

IV.  SURVEY  DATA  ANALYSIS 12 

V.  STATISTICAL  MODEL  TO  PREDICT  DATA  GAPS 2 5 

VI.  EVALUATION  OF  SURVEY  RESULTS 4 8 

VII.  CONCLUSION 51 

APPENDIX  1 - SURVEYOR  DATA 1-1 

APPENDIX  2 - MBTA  FILE  STRUCTURE 2-1 

APPENDIX  3 - PROGRAM  EXTRAC 3-1 

APPENDIX  4 - PROGRAMS  EDIT  AND  EDITl 4-1 

APPENDIX  5 - PROGRAMS  EDIT2  AND  EDIT3 5-1 

APPENDIX  6 - PROGRAM  SEPAR 6-1 

APPENDIX  7 - FINAL  OUTPUTS:  PAKRIV  AND  PAKLEC..  7-1 

APPENDIX  8 - PROGRAM  PROPFL 8-1 

APPENDIX  9 - STATION  PROFILE:  PAKPFL 9-1 

APPENDIX  10  - PROGRAM  CONTRAC 10-1 

APPENDIX  11  - PROGRAM  COMBIN 11-1 

APPENDIX  12  - PROGRAM  ESTRAG 12-1 

APPENDIX  13  - ANOVA  RESULTS 13-1 


V 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1.  SCHEDULING  OF  THE  PASSENGER  SURVEY 5 

2.  DISTRIBUTION  OF  LOADING  PASSENGERS:  CASE  1 17 

3.  DISTRIBUTION  OF  UNLOADING  PASSENGERS:  CASE  1 18 

4.  DISTRIBUTION  OF  LOADING  PASSENGERS:  CASE  2 2 0 

5.  DISTRIBUTION  OF  UNLOADING  PASSENGERS:  CASE  2 21 

6.  DISTRIBUTION  OF  LOADING  PASSENGERS:  CASE  3 2 3 

7.  DISTRIBUTION  OF  UNLOADING  PASSENGERS:  CASE  3 24 

8.  PLOT  OF^THE  RESIDUALS,  Eij,  VS.  THE  PREDICTED 

VALUE,  Rij 34 

9.  NORMALITY  PLOT  OF  THE  STANDARDIZED  RESIDUALS 

(UNTRANSFORMED  DATA)  36 

10.  PLOT  OF  PREDICTED  VALUES,  Rjj,  VS.  THE  ACTUAL,  Rt y 

(UNTRANSFORMED  DATA) T . . 37 

11.  PLOT  OF  THE  RESIDUALS,  Ejj,  VS.  THE  PREDICTED 

VALUES,  Rij  (TRANSFORMED  DATA) 45 

12.  NORMALITY  PLOT  OF  THE  STANDARDIZED  RESIDUALS,  Ejj 

(TRANSFORMED  DATA) 46 

13.  PLOT  OF  THE  PREDICTED  VALUES,  R..,  VS.  THE  ACTUAL, 

Rij 47 

14.  PROFILE  OF  LOADING  PASSENGERS  AT  PARK  STREET 

STATION,  1978  (TO  RIVERSIDE) 52 


VI 


LIST  OF  TABLES 


Table  Page 

1.  SURVEY  PLAN 7 

2 . MBTA  FILE  STRUCTURE 9 

3.  EDITING  PROGRAMS 10 

4.  OBSERVED  PASSENGER  LOADING  RATES  AT  STATIONS 

(TRAINS  TO  RIVERSIDE) 13 

5.  OBSERVED  PASSENGER  UNLOADING  RATES  AT  STATIONS 

(TRAINS  TO  RIVERSIDE) 13 

6.  OBSERVED  PASSENGER  LOADING  RATES  AT  STATIONS 

(TRAINS  TO  NORTH  STATION) 14 

7.  OBSERVED  PASSENGER  UNLOADING  RATES  AT  STATIONS 

(TRAINS  TO  NORTH  STATION) 14 

8.  ANOVA  RESULTS  ON  RIVON.MTY;  MODEL  (2)  ANALYSIS  OF 

VARIANCE  FOR  THE  REGRESSION  (UNTRANSFORMED  DATA).  32 

9.  ANOVA  RESULTS  ON  RIVON.MTY:  MODEL  (3)  ANALYSIS  OF 

VARIANCE  FOR  THE  REGRESSION  (TRANSFORMED  DATA)...  39 

10.  ESTIMATED  LOADING  RATE  (ALL  RIVERSIDE  TRAINS) 49 

11.  ESTIMATED  UNLOADING  RATES  (ALL  RIVERSIDE  TRAINS).  49 

12.  ESTIMATED  LOADING  RATES  (ALL  NORTH  STATION 

TRAINS) 50 

13.  ESTIMATED  UNLOADING-  RATES  (ALL  NORTH  STATION 

TRAINS) 50 


vii/viii 


I.  Introduction 


This  paper  presents  a sampling  scheme  and  an  estimation 
methodology  to  approximate  passenger  flows  at  a system  of 
underground  stations.  The  resulting  data  will  be  calibrated 
into  a simulation  model  whose  objective  is  to  optimize  the 
scheduling  of  the  Light  Rail  Vehicle  (LRV)  to  meet  the  daily 
changing  passenger  demands.  Such  an  effort  is  part  of  an 
operation  analysis  system  designed  by  the  Transit  Systems  Branch 
Urban  Systems  Division,  Office  of  Ground  Systems,  to  assist 
transit  properties  in  their  management  of  transit  systems 
maintenance  and  operations. 

When  the  LRV's  were  introduced  in  January  of  1977  to 
Boston’s  Massachusetts  Bay  Transit  Authority  (MBTA)  to 
gradually  replace  the  old  PCC’s  (President  Conference  Com- 
mittee) on  the  Green  Line,  it  was  deemed  necessary  by  the 
U.S.  Department  of  Transportation,  who  funded  80%  of  the 
Boeing  made  cars,  to  identify  areas  of  operational  improve- 
ment and  to  conduct  cost  reduction  studies.  The  project,  spon- 
sored by  the  Department’s  Urban  Mass  Transportation  Administra- 
tion, aims  at  employing  simple  analytical  and  data  management 
tools  to  efficiently  develop  a prototype  operations  analysis 
system  which  can  be  modified  and  applied  to  various  transit 
properties  to  fit  their  specific  transit  management  needs. 

Specifically,  using  MBTA’s  Riverside  Line  as  the  experi- 
mental unit,  this  study  collects  and  analyzes  passenger  data 
at  the  thirteen  underground  stations.  Thirteen  daily  profiles 
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of  passenger  loading  and  unloading  rates  are  estimated  so  that 
the  distribution  of  passenger  flow  at  any  given  station  at  any 
time  can  be  obtained.  By  varying  the  headways  (and  thus  the 
schedules)  of  the  trains,  one  can  arrive  at  the  optimal  train 
schedule  to  accommodate  the  expected  passenger  demand. 

II.  The  Survey  Plan 

A.  Knowledge  gained  from  the  Survey  at  the  Surface  Stations 

In  May  1977,  a sample  of  LRV  trips  to  and  from  Riverside 
was  selected  and  the  loading  and  unloading  passenger  counts 
were  taken  by  an  onboard  surveyor.  This  was  done  only  for 
that  part  of  the  trip  when  the  train  remained  on  the  surface 
and  when  the  patrons  of  the  LRV  were  allowed  to  board  and 
alight  through  only  one  door.  The  underground  stations  pose  a 
more  complex  problem.  The  am.ount  of  patronage  increases  several 
folds;  all  the  trains  are  merged  onto  one  track;  and  at  least 
three  doors  are  open  as  a train  pulls  in  the  station.  A more 
efficient  counting  procedure  is  necessary,  so  that  the  distri- 
bution of  passenger  flow  at  these  stations  for  the  Riverside 
trains  can  be  estimated  with  some  degree  of  confidence  without 
the  investment  of  a large  amount  of  resource  and  time. 

Results  of  the  analysis  of  1977  surface  station  data  have 
proved  useful  not  only  in  the  delineation  of  the  sampling  plan 
for  the  underground  stations,  but  also  in  the  verification  of 
the  resulted  passenger  profile  thus  obtained.  First,  one 
major  finding  in  the  analysis  indicates  an  insignificant  day 
to  day  difference  in  the  total  passenger  volume  using  the  LRV. 


Second,  it  was  also  confirmed  by  statistical  chi-square  tests 
that  each  surface  station  holds  in  general  a constant  market 
share  of  the  total  passenger  volume  both  for  the  morning  and 
the  afternoon  hours.  Hence  it  was  possible,  given  a total 
trip  passenger  volume  at  any  time,  to  ascertain  an  expected 
passenger  volume  from  a certain  station  on  the  way.  It  is 
safe  to  deduce,  therefore,  that  "station’’'  (or  location)  and 
"time  of  day"  are  the  two  most  important  variables  which  deter- 
mine potential  passenger  flows  onto  the  LRV’s.  That  is,  we  pro- 
pose to  investigate  how  the  rate  of  passenger  flow  varies 
systematically  with  the  two  factors.  The  length  of  the  head- 
way, another  essential  parameter,  will  of  course  determine  the 
actual  congestion  and  the  waiting  time. 

B.  The  Survey  Plan  at  the  Underground  Stations 

In  May  1978,  a new  survey  plan  was  devised  whereby 
H.  H.  Aerospace,  Inc.,  a minority  contractor,  would  send  a team 
of  eight  surveyors  to  the  thirteen  underground  stations.  Their 
mission:  to  collect  passenger  data  on  the  Riverside  trains, 

both  inbound  and  outbound. 

It  is  necessary  to  discuss  the  survey  plan  and  field 
assignment  in  some  detail  because  it  is  precisely  the  design 
of  the  survey  that  enables  the  analyst  to  draw  meaningful 
inference  from  the  survey  data.  The  objective  of  the  survey 
is  to  derive  a representative  daily  distribution  of  passenger 
flows  either  by  summarizing  the  survey  data  or  by  means  of 
estimation  using  a simple  statistical  model.  The  use  of  a 
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statistical  model  economizes  the  number  of  observations  re- 
quired, because  variation  in  the  data  is  minimized  by  the 
control  variables  in  the  model.  For  example,  indicator  varia- 
bles such  as  the  location  of  the  station  and  the  time  of  the 
day  probably  explain  much  of  the  data  variance,  and  a two  factor 
analysis  of  variance  (ANOVA)  model  seems  most  appropriate.  The 
implementation  of  the  model  is  the  major  criterion  for  the 
survey  design. 

Two  other  constraints  are  the  time  span  (five  days)  in 
which  the  survey  is  to  be  conducted,  and  the  labor  resource 
(eight  surveyors).  The  number  of  surveyors  assigned  to  a 
station  depends  on  the  amount  of  passenger  traffic  the  station 
receives.  For  example,  four  surveyors  are  assigned  to  Park 
Street  Station  at  peak  hours,  while  Auditorium  or  Boylston  has 
been  adequately  manned  by  only  one  surveyor.  Figure  1 dis- 
plays a matrix  of  110  cells,  40  of  which  are  randomly  selected 
and  assigned  to  a team  of  surveyors  (numbering  from  1 to  4 
depending  on  the  needs  of  the  station/shift.)  There  are  two 
shifts  for  each  station  and  each  direction.  Science  Park  and 
Lechmere,  being  the  stations  near  or  at  the  end  of  the  line, 
have  surveyors  only  at  the  outgoing  (towards  Riverside)  direc- 
tion. The  random  selection  of  the  working  shifts  forms  the 
primary  sampling  stage.  The  complexity  of  the  plan  arises 
when  only  eight  people  are  sent  to  the  field  to  work  in  these 
shifts.  Here  the  objectivity  of  the  survey  design  gives  way 
to  the  subjective  selection  of  the  time  span  within  each  shift 
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FIGURE  1.  SCHEDULING  OF  THE  PASSENGER  SURVEY 


at  various 


based  on  a priori  knowledge  of  the  passenger  flow 
stations  and  the  desire  to  sample  more  heavily  the  busy  sta- 
tions and  the  peak  hours.  Each  station  is  to  be  observed  at  one 
time  or  another  and  no  time  period  is  missed.  The  result  is  the 
survey  plan  delineated  in  Table  1,  which  gives  Kenmore  contin- 
ual coverage  as  a point  of  reference  to  other  stations. 

The  arrangement  of  labor  resource,  the  suggestion  of  an 
optimal  position  from  which  the  sampler  can  have  a good  view 
of  the  car  doors,  the  layout  of  the  station  platform,  and 
specific  instructions  during  peak  hours  are  all  part  of  a two 
hour  or ientat ion/ training  for  the  samplers.  Where  it  is 
difficult  to  count  passengers  at  all  doors  (for  example, 
during  rush  hours  at  busy  stations)  the  samplers  are  to 
randomly  select  a car  door  for  the  count  but  systematically 
go  to  the  next  door  of  the  second  train  when  it  comes  around. 
This  ensures  the  representativeness  of  the  data  when  the  counts 
are  extrapolated  to  represent  the  entire  train. 

Appendix  1 is  a sample  of  the  data  reported  by  a surveyor. 
Specifically,  the  data  elements  are  the  stations,  the  time  of 
train  arrival,  the  direction  of  the  trip,  the  number  of 
loading  and  unloading  passengers  at  the  car  door,  and  the  car 
characteristics.  The  samplers  are  also  to  report  on  other 
trains  that  pull  in  the  station  without  trying  to  measure 
their  loading  and  unloading  activities,  and  note  any  un- 
usual circumstances  such  as  large  group,  train  malfunction, 
etc . 
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TABLE  1.  SURVEY  PLAN 


of 

Station 

Monday 

Tuesday 

Wednesday 

Thursday 

Friday 

AM 

Kenmore 

PM 

X 

0 

X 

0 

X 

AM 

Auditorium 

PM 

X ® 

X 

0 

AM 

Copley 

PM 

X 

X 

AM 

Arlington 

PM 

(D 

X 

X 

AM 

Boylston 

PM 

X 

0 

X ® 

AM 

Park  St. 

PM 

X 

X 

AM 

Gov't  Center 

PM 

X 

X 

0 0 

AM 

Haymarket 

PM 

X 

X 

0 © 

X 

AM 

North  Station 

PM 

X 

X ® 

AM 

Science  Park 

PM 

AM 

Lechmere 

PM 

X 

X Riverside  to  Lechmere  (N.S.) 
® Lechmere  to  Riverside 
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To  confirm  the  notion  of  nonsignificant  daily  variation 
of  passenger  demand,  the  survey  design  incorporates  opportun- 
ities to  retest  that  hypothesis.  Samples  are  taken  during  the 
same  time  span  for  different  days  of  the  week  at  Auditorium, 
Arlington,  and  Park  Street. 

The  design  in  effect  offers  the  analyst  an  unreplicated 
discontinued  time  series  of  passenger  flow  for  each  station, 
unreplicated  because  no  time  periods  (except  those  for  hypothe- 
sis testing  purposes)  are  repeated  more  than  once,  and  dis- 
continued in  the  sense  that  data  gaps  exist  due  to  the  labor 
constraint.  It  is  the  intention  of  Section  IV,  however,  to 
attempt  to  fill  these  data  gaps  with  an  adequate  statistical 
model  and  to  estimate  the  daily  trend  of  passenger  demand 
with  some  degree  of  statistical  confidence. 

III.  Data  Editing 

Like  any  other  data  collection  activity,  editing  and  data 
screening  becomes  an  essential  part  of  the  total  effort  to 
prepare  the  data  for  analyses.  Anomalies  or  obvious  errors 
are  to  be  corrected  or  discarded.  If  more  than  one  sur- 
veyor is  surveying  the  train,  each  being  responsible  for  a car 
or  a door,  the  results  of  each  have  to  be  combined.  Or,  if  the 
combined  passenger  count  still  does  not  cover  the  entire  train, 
extrapolation  procedures  are  to  be  developed. 
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When  the  survey  is  completed,  the  data  sheets  are  de- 
livered to  the  keypunchers  for  the  transcription  of  data 
into  a file  area  in  DEC-10.  The  file  structure  is  outlined 
in  Table  2. 


TABLE  2.  MBTA  FILE  STRUCTURE 


Field  Name 
MAYDAY 
TIME 
LINE 


TRAIN  SIZE 
TYPE  OF  CAR 


PASSON 

PASSOFF 


Field  Description 
Survey  date 
Time  of  train  arrival 
Line  of  train: 

R=Riverside 
CC=Cleveland  Circle 
BC=Boston  College 
H=Hunt ing ton/ Arborway 
K=Kenmore 
N=North  Station 
L=Lechmere 

GC=Government  Center 

P=Park  Street 

Number  of  cars  in  the 
train 

LRV=Light  Rail  Vehicle 

PCC=Pres ident ' s Conference 
Committee 

Number  of  loading  pas- 
sengers 

Number  of  unloading 
passengers 


Data  Type 
Date 
Integer 
Text 


Integer 

Text 


Integer 

Integer 


Length 

8 

4 

2 


1 

3 


3 

3 
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TABLE  2 (CONT.) 


Field  Name 

Field  Description 

Data  Type 

Length 

CARl 

X if  car  number  1 only 
or  if  the  trainsize  = 1, 
blank  if  otherwise 

Text 

1 

CAR2 

X if  car  number  2 is 
counted 

Text 

1 

CAR3 

X if  car  number  3 is 
counted 

Text 

1 

DOORl 

X if  door  is  counted, 
blank  otherwise 

Text 

1 

D00R2 

X if  door  2 is  counted, 
blank  otherwise 

Text 

1 

D00R3 

X if  door  3 is  counted, 
blank  otherwise 

Text 

1 

COLLECT 

Surveyor's  identifier 
from  1 to  4 since  a 
maximum  number  of  4 
workers  can  be 
assigned  to  a shift 

Integer 

1 

In  the  file,  each  entry,  representing  one  recording  made  by 


one  surveyor  has  13  fields.  A series  of  editing  programs  delineated 


in  Table  3 brings  the  file  to  its  final  form,  with  each  entry  repre- 


senting one  train  observation  at  a station.  (See  Appendix  7.) 

TABLE  3.  EDITING  PROGRAMS 


Program  Name  Function 


Output  File 


EXTRAC 
(Appendix  3) 


Extracts  individual  station  station  file 
data  from  the  main  MBTA  file;  such  as 
sorts  data  by  date  and  time,  "ARLINGTON"  or 
line,  etc.,  so  that  the  "PARK  ST" 

working  area  is  reduced. 

It  is  here  where  data  are 
screened  manually  for 
obvious  errors. 


EDIT  Uniquely  identifies  each 

(Appendix  4)  car  door  whose  passenger 

count  is  included  in  the 

surveyor's  record. 


STRAl  (working 
file) 
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TABLE  3.  (CONT.) 


Program  Name 
EDIT2 

(Appendix  5) 


SEPAR 

(Appendix  6) 


PROFL  , . 
(Appendix 


8) 


CONTRA 

(Appendix  10) 


Function 

Takes  STRAl  as  input , 
combines  passenger  counts 
from  each  surveyor,  extra- 
polating if  necessary,  to 
derive  total  train  count* 

Separates  STRA2.DMS  into 
2 files  according  to  the 
direction  of  the  train 
(inbound  or  outbound)  and 
within  each  file  sorts  the 
data  by  time. 

Takes  PAKLEC  or  PAKRIV  and 
computes  headways,  loading 
and  unloading  rate.  Also 
counts  the  number  of  trains 
between  Riverside  trains. 


Takes  ARLPFL  or  ARLPFL. 
NOR  and  computes  loading 
unloading  rates  for  each 
hour  interval. 


Output  File 

STRA2 .DMI 
(working  file) 
to  be  made  into 
a 1022  data  file. 
(STRA2.DMS) 

eg.  PAKLEC  and 
PAKRIV  (final 
data,  two  files 
for  each  station) 
(Appendix  7) 

eg.  PAKPFL 
(Riverside)  and 
PAKPFL. NOR  (North 
Station)  or  PAKPFL. LEC 
(Lechemere) . Eacn 
file  represents  an 
average  daily  pas- 
senger profile  (raw 
data)  at  the  station. 
(Appendix  9) 

PAKCON.RIV  or 
PAKCON.NOR  (working 
files) 


and 

half 


COMBIN 

(Appendix  11) 


Combines  all  ***CON.RIV  or  RIVON.MTY 

***CON.NOR  files  into  two  RIVOFF.  MTY 

data  matrices  (one  for  load-  or  NORON.  MTY 
ing  passenger  rates  and  one  NOROFF . MTY 
for  unloading  passenger  (Table  2) 

rates)  with  "Stations"  as 
the  Columns  and  "time 
periods"  as  the  Rows 


*The  extrapolating 
TOTAL  Train  Count  = 


method  is; 

total  passengers  count 
number  of  doors  counted 


X 


total  no.  of  doors 
on  a train 


Assumes  that  passengers  board  trains  equally  at  all  doors  which 
is  not  true  for  certain  stations  during  certain  time  perio  s. 
Passengers  tend  to  board  door  nearest  turnstile  entrance  or 
stairwell . 
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The  last  two  programs  are  necessary  to  prepare  the  data 
for  input  into  a statistical  model  which  synthesizes  the  in- 
formation for  all  stations  and  all  time  periods  to  predict 
missing  values,  derives  general  trend  for  each  station  and 
estimates  the  average  variation  of  passenger  demand.  Such  a 
model  is  discussed  in  the  next  section.  The  matrices  in  Tables 
4 through  7 show  not  only  what  is  missing  but  give  a bird's- 
eye-view  of  the  sampling  design. 

IV.  Survey  Data  Analysis 

One  of  the  assumptions  underlying  the  sampling  design  is 
that  one  can  pool  the  observations  made  at  different  days  of 
the  week  to  form  a representative  daily  profile  of  an  average 
week  day.  This  assumption  eliminates  the  necessity  of  obtain- 
ing sample  data  for  each  day  of  the  week  and  deriving  five 
daily  profiles  for  each  station.  Previous  analyses  on  the 
surface  stations  reflect  little  day  to  day  difference  in  total 
passenger  traffic.  Our  sampling  design  for  the  underground 
stations  also  allows  for  similar  testing  of  the  hypothesis. 

This  section  presents  three  cases  in  which  two  independent 
sets  of  passenger  counts  are  obtained  at  a station  during 
the  same  time  span  but  for  different  days  of  the  week.  To 
each  case,  a nonparametr ic  test,  namely,  the  Mann-Whitney 
U test  is  applied  to  see  if  the  distribution  of  the  observa- 
tions taken  on  Day  1 is  different  from  those  taken  on  Day  2. 
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TIME 

700 

7 3'J 

Boo 

630 

930 

lOOO 

1 030 

1100 

1130 

1200 

1230 

1300 

l33o 

1400 

1430 

15  00 

l5  3ij 

l60o 

lb30 

1700 

1730 

1800 

l83o 

l900 

1930 

2000 

2030 


TABLE  4.  OBSERVED  PASSENGER  LOADING 

RATES  AT  STATIONS  (TRAINS  TO 
RIVERSIDE) 


0,00 

0,00 


0,00 

0,00 


0,50 
n . n (■» 


0,00 


0,00  0,00 
0,00  ■ 0,00 


0,00 

0,00 


HAY. 

N. 

STA. 

0,00 

0,00 

2,24 

0,00 

4,87 

ITTVS 

2,06 

13,19 

5,22 

3,90 

0,00 

4,60 

T7TT 

3,90 

2,22 

O,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0.00 

0,00 

2,94 

0,0  0 

1,00 

1,00 

1,75 

1.93 

2,77 

2,45 

0,00 

l,b9 

0,00 

3,84 

0,00 

1,65 

0,00 

2,15 

0 , 00 

0,00 

0,00 

EH] 

0,00 

. 0,00 

0,00 

: U.UO 

0,00 

: 0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

SCI. 

LECH. 

PK. 

0,00 

0,00 

0,00 

0,00 

0,30 

0,0  0 

0,28 

0,00 

0,21 

0,00 

0.16 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0.00 

0,00 

i.7l 

0,00 

3,58 

0,00 

6,16 

0,00 

3.47 

0,00 

4.52 

0,00 

5.70 

0,00 

0,00 

0,00 

0.00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0.00 

0,00 

0,00 

0,00 

0,00 

TABLE  5. 


OBSERVED  PASSENGER  UNLOADING  RATES 
AT  STATIONS  (TRAINS  TO  RIVERSIDE) 


TIME 

KEN. 

AUD. 

COP. 

ARL. 

BOYL. 

PARK 

GOV. 

HAY. 

1 

0,00 

CTR. 

700 

0,00 

0,00 

4,57 

0,00 

0,00 

0,00 

0,00 

730 

0,0  0 

o.ooi 

3.16 

o.oo! 

0,00 

1.22 

0.00 

0,U4 

80o 

3.00 

0,00 

6,b8 

0,00 

1.97 

TT^ 

0,09 

830 

4,19 

0,00 

8.19 

7.91 

0,00 

4,74 

1,96 

0,06' 

900 

5.05 

0,00 

0.00 

4.26 

0,00 

3,20 

0,57 

0,13 

930 

1.79 

0,00 

0,00 

2.17 

0.00 

1,76 

0.44 

0,00 

1 OOO 

1.33 

T73T 

0,00 

0,00 

0,61 

3.25 

0,81 

0,00 

10  30 

1,44 

0,b6 

0,00 

77^ 

1.65 

1,10 

0,43 

0,00 

1100 

1.78 

0.35 

0,00 

0,55 

0,20 

0.00 

0,00 

0,00 

1 130 

1.55 

0,71 

0,00 

0,60 

0.89 

0.00 

0,00 

o.oo' 

1200 

1 ,b8 

0,73 

0.00 

2.63 

0,0  {T 

iO.OO 

\ 

0,00 

o.uo! 

1230 

C.48 

0.00 

6.00 

1^51 

0,00 

10,00 

0,00 

0,00 

1300 

0,00 

0,00 

' 

3.03 

0.00, 

0,00 

10,00 

0,00 

0,00 

l330 

1.7T 

0,00 

2,80 

O.OOj 

0.00 

0.00 

0,00 

0.00 

1400 

2.32 

0,00 

' 

3.15 

0,00; 

77^ 

ITT? 

0,00 

77T^ 

1430 

6,42 

10,00 

^.20 

0,00 

0.18 

2.28 

' 

0,00 

0,07 

1500 

3.81 

'o,oo 

0,00 

0,00 

0,34 

1.87 

0,00 

0,03 

1530 

1,87 

0,00- 

0,00 

0,00 

0,20 

7.00 

0,00 

0.38 

1600 

2,66 

0,00 

0,00 

0,00 

2.60 

3.38 

77T7 

0,00 

163o 

3.38 

0,00 

0.00 

0,00 

1 ,09 

4.31 

1.82 

OtOQ 

1700 

7,38 

0,00 

0,00 

0,00 

1.34 

8,76 

1,05 

1730 

10,04 

0,00 

0,001 

0,00 

U.18 

5.37 

0,87 

0.00 

1800 

6,o6 

0,00 

0.00 

0.00 

0^54 

4,00 

0,63 

0,50 

1 830 

2,75 

0,00 

0,00 

0.00 

5,70 

0,44 

o,o6( 

1900 

o,0o 

0,00 

o.oo 

0,00 

0,00 

o76fi 

0,07 

0,00 

1930 

0,00 

0,00 

0,00 

0.00 

0.00 

0,00 

0,00 

0,00 

2000 

0,00 

0,00 

1 .50 

0,00 

0,00 

0.00; 

0,00 

o.oo! 

203o 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

'o.oo 

o.od] 
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TABLE  6.  OBSERVED  PASSENGER  LOADING  RATES 

AT  STATIONS  (TRAINS  TO  NORTH  STATION) 


TIME 

KEN. 

AUD. 

COP. 

ARE. 

BOYL. 

PARK 

GOV. 

HAY. 

1 

i 

CTR. 

700 

0 ,00 

'0,00 

0,33 

0,00 

0,00 

0.00 

10,00 

0.00; 

730 

0,00 

:0,00 

0,47 

jO.OO 

o.oo' 

0,00 

0,00. 

O.OOj 

800 

o,00 

jo, 00 

1.75 

0,00 

0,00 

T7T0 

o;i4 

0,00, 

830 

3,06 

10,00 

3.36 

1,84 

0,00  1 

2.28 

0,13 

O.OOj 

900 

5,01 

!o,oo. 

0,00 

2.17 

0,00  , 

0.63 

0,43 

0,00 

930 

1.14 

o.oo; 

0.00 

0.42 

0,00  , 

0,79 

0,54 

77n 

lOOO 

0,35 

0,00 

0,00 

0,33 

1.21 

0,15 

0,00 

!030 

1.14 

1 

1,66 

0,00 

0,92 

0,07 

1,00 

0,33 

0,Q6 

1100 

3,29 

1.13 

OTT? 

1,00 

0,16 

0,81 

0.94 

0,00 

1 130 

0,93 

1.45 

1.26 

0,36 

0,62 

0,00, 

0,11 

0,00 

1 200 

1.75 

0.00 

0,46 

0,56 

bTW 

0,00, 

0,10 

0,00 

1230 

1 ,66 

O.OOj 

5.23 

0,00: 

0,00 

0,00; 

0,03 

0,00 

1 300 

1.11 

0.00; 

0 , 0 b 

O.OOi 

0,00 

0,00, 

0,00 

1330 

2,94 

0,00 

O.OOi 

0,00 

O.OOi 

0,00 

o.oT 

1400 

3j84 

9. 00 

2j38 

9 .00 

[0,00 

o,oo( 

0,00 

0,04 

1430 

3,96 

0,00 

1.'74 

0,00 

^0,00 

0.00; 

0,00 

0,05 

1500 

3,00 

9.00 

3,11 

0,00 

0,00 

0.00 

0,00 

6,  Ob' 

1530 

0,00 

0.00 

0,00 

10.00 

0,00 

5.00 

0,00 

0,00 

1600 

2,00 

9,0^. 

0,00 

0,00,  1 

0.00 

1,80 

1737 

0.00 

1630 

3,67 

2.54 

0,00 

0,00 

Fin 

5,81 

0,24 

FTIT 

1700 

; 

1.36 

0,00 

1,58 

1,09 

.2.18, 

0j90 

2.22 

1730 

4,35 

4,90 

o.cb 

3.90 

0.53 

5,70 

0,97 

0.24 

1800 

1,08 

1.42 

0,00 

2,50 

0,08 

i.oo. 

0.72 

0.00, 

l830 

0,03' 

:0,00 

lo.oo 

1 

1.1? 

0.25 

liJl 

0,26 

O.OOj 

1900 

^00 

0,00 

10,00; 

.0,00 

0,00 

0,00. 

FTCnr 

0,00. 

1930 

0,00 

0.00 

|0,00 

O.OOi 

0,00 

0,00 

0,00 

|0,00 

200o_  io,oo 

0,00. 

.OjOO 

!o,oo 

lOjOO 

QjOO. 

0.00 

[0^0^ 

2030 

0,00j 

0.00 

0.00' 

0,00 

'o.oo 

0.00 

'b.bo 

(0,00 

TABLE  7.  OBSERVED  PASSENGER  UNLO'ADING 
RATES  AT  STATIONS  (TRAINS  TO 
NORTH  STATION) 


TIME 

KEN. 

AUD. 

COP. 

ARE. 

( 

700 

v^'.OO 

0,00 

0,47 

[0,00 

7 3n 

J,0  0 

0 ,0V 

0.95 

'o.oo! 

800 

. 0 0 j 

0,0  0 
o.oo' 

2,69 

0,00 

830 

w ,50 

3.96 

3.90 

900 

v,30 

1 

0,00 

O.OOi 

4,67 

930 

! .59 

0,00 

C 

• 

c 

c 

2,9b 

1000 

0.35 

llibT 

O.OOi 

0,00 

1030 

. .47 

1.45 

0,00' 

1 ,46 

11  Oo 

V.  .43 

1.45 

1.44 

1.57 

1 130 

J .00 

y , 00 

0,85 

0,62 

1200 

^ .32 

0,00 

1.19 

U.96 

l23y 

V.  .88 

0.00 

1.15 

0,00; 

1300 

.33 

io » 0 0 

0,00 

O',  00  i 

1330 

.50 

lo.oo 

0,68 

O.OOi 

1 400 

1 ,o9 

lo.oo 

2.00 

0 , 0 0 i 

14  30 

1.04 

10.00 

0,72 

O.OOi 

1500 

•J  . C 0 

lo.OO ; 

0 . 86 

o.oo; 

1530 

V . OO 

lo.oo 

; O.UC- 

0,00: 

1600 

>.29 

0 * 0 C' 

>0,00 

; 0.00 

16  30 

-.3  9 

oT?5' 

■0,00 

0,00 

1700 

-.23 

0,71 

O.oo 

0,50 

17  30 

:/,39 

0,90 

0,00 

0,94 

l800 

■;,5  0 

0.77 

0,0.0 

0,8  8 

1830 

1,37 

0,00 

0,00 

0.05 

1900 

■ ; ,00 

0,00 

'0,00 

0 ,oO 

1930 

'«■ , 0 0 

0,00 

0,0  0 

0,00 

2000 

■; . 0 0 

iO  , OO 

0,0  0 

0,00 

2030 

'J  . 0 0 

lo.oo 

o 

c 

0.00 

BOYE. 

PARK 

GOV. 

HAY. 

N. 

CTR., 

STA.  1 

0,00 

0,00 

0,00 

0,00 

o,oo| 

0,00 

o 

O 

C 

' o.oo' 

0,00 

io.oo’ 

0,00' 

22, U7 

20,43 

0,00' 

5,15 

0,00* 

10,69 

11,65 

o.oo! 

1.09 

0,00 

13.07 

16,97 

0,00 

3.53 

0,00' 

5,56 

7,15 

0,47 

1,09 

brn 

8,66 

6,88 

0,94 

o,6o' 

0.52 

5.11 

2,93 

2,22 

'0,00| 

O.bl 

7.19 

3,97 

,0.49 

O.OOi 

1.12 

O.oo; 

1.14 

0,00; 

0,00: 

0 ,00_ 

O.OOi 

5,10 

o.ooj 

0,00; 

0,00 

0 , 0 0 j 

4,92 

'0,00' 

0,00: 

0,00 

o,uo; 

0,00 

0,00. . 

3,56 

0,00' 

o.oo! 

; 0,00 

0,50 

2.14 

0,00 

o.oo' 

0,00 

2.39 

1,82 

0,00 

■ o.ooj 

! 0,00, 

2,26 

2,36 

0,00 

0.00 

0,00 

3.69 

0,00 

0,00 

10,33 

0,00 

0,00 

'O.ool 

0,00 

8,17 

5.52 

0,00 

10,00 

oTW 

16.22 

4.79 

3.96 

0.00, 

0.74 

l0,30 

4,67 

9,04 

10, 00' 

0.41 

12.33 

5,34 

9,05 

•0,00 

0,46 

7,76 

3.62 

3,41 

0,00 

0,25 

3.66 

2,39 

0,00 

0,00 

0,00 

0.00 

0,00 

i 0,00 

0,00 

0,00 

0,00 

0,00' 

0,00 

0,00 

0,00 

O.OO- 

0,00 

0,00 

•0.09 

0,00 

0,00 

0,00 

0,00 

0,00 
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The  method  requires  the  assignment  o£  rank  to  each 
of  the  n^  + n^  observations  according  to  its  magnitude,  where 
n^  and  n^  are  the  total  numbers  of  observations  for  Day  1 
and  Day  2 respectively.  A U statistic  is  then  formed  by 

. . . "A"! 

''l"2  2 

i €:  Day  1 


U'  = "i"2  " 


- 1) 


Z 


W. 


i t Day  2 


A two  tailed  test  for  the  hypothesis,  H,  that  the  two 
samples  come  from  the  same  population  and  thus  have  the  sam.e 
distribution  is  as  follows: 

At  a significance  level  a = .05,  is  rejected  if 


min  (U,U')  < ,np 

min  (U,U')  > 


or 


where  C 


(a/2,  n^,n2)  is  obtained  from  the  table  of  critical 
values  of  U in  the  Mann-V/hitney  Test  (Appendix  14)  and 


‘'(1-0/2, n^n^)  " "l"2  ' ^ (o/ 2 ,n  ^ .n^)  • 


In  all  of  the  three  cases  in  the  following,  the 
test  confirms  the  assumption  of  insignificant  day  to  day 
variation  in  the  average  passenger  flow  for  the  underground 
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stations.  Along  with  each  test  result,  the  actual  data  is 
plotted  graphically  to  reflect  the  similarity  of  their 
d i St r ibut ions . 

Case  1:  Park  Street  Station,  l;30pm  - 3:30  pm  (Figures  2 and  3). 

Loading  Passengers 


U = 16  X 16  + 


16(16  + 1) 
2 


224 


168 


U’  = 16  X 16  - 224  = 88 


min  (U,U')  = 88 


C 


a/  2 


75 


Hence  ^ min(U,U’)  < ^o  accepted. 


Unloading  Passengers 


U = 16  X 16  + 


16(16  + 1) 
2 


279.50  = 112.50 


U*  = 16  X 16  - 112.50  = 143.50 


min  (U,U')  = 112.50 


C 


(a/2,n^,n2)  = 75 


*^(l-a/2  ,n^  ,n2)  = 181 


Hence  C 


and  H^  is  accepted. 
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FRIDAY 


THURSDAY 


Passengers  on  Passengers  on 

n2=16  n2=16 


• Trains  other  than  Riverside  Trains 


FIGURE  2.  DISTRIBUTION  OF  LOADING  PASSENGERS:  CASE  1 
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FRIDAY 


THURSDAY 


1 : 30  pm  1 : 30  pm 


No.  o£  Passengers  Unloading  No.  Passengers  unloading 


112  = 16 


n^  = 16 


» Trains  other  than  Riverside  Trains 


FIGURE  3.  DISTRIBUTION  OF  UNLOADING  PASSENGERS:  CASE  1 
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Case  2:  Auditorium  Station,  10:00  am  - 12:00  pm  (Figures  4 and  5). 

Loading  Passengers 

U = 10  X 13  + - 116.5 

= 68.50 

U’  = 61.50 

min  (U,U’)  = 61.50 

->  = 33 
Ca/2 ,n^ ,n^) 

Cri  /o  . = 97 

(1 -a/2 ,n^ 

Hence  ,n^  1 rain  (U,U')  < C ^ ’"l ’"2^ 

and  is  accepted. 

Unloading  Passengers 

U = 10  X 13  + - 118 

= 67 
U’  = 63 

min  (U,U’)  = 63 
^(a/2,n^,n2)  " 

^(l-a/2 ,n^ ,n2)  = 97 


Hence  C 


(1 -a ,n^ ,n^) 


£ min  (U  ,U ’ ) < C 


(1 -a/2 ,n^ 


FRIDAY 


TUESDAY 


n2=13  n^=lC 


• Trains  other  than  Riverside  Trains 


FIGURE  4.  DISTRIBUTION  OF  LOADING  PASSENGERS:  CASE  2 
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FRIDAY  TUESDAY 


* Trains  other  than  Riverside  Trains 


FIGURE  5.  DISTRIBUTION  OF  UNLOADING  PASSENGERS:  CASE  2 
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Case  3:  Arlington  Station,  10:00  am  - 12:00  pm  (Figures  6 and  7) 


U 


8 X 11  + 


8(8  + 1) 
2 


50 


74 


U'  = 8 X 11  - 50  = 38 

min  (U,U')  = 38 

->=19 
(a/2 ,n^ ,n^) 


C r,  / T -.  = 69 

(l"a/2  ,n^^  fTL 2^ 


Hence  C 


(a/2,n^,np  i < 


and  is  accepted 
Unloading  Passengers 


U = 8 X 11  + - 89 

= 35 

U’  = 8 X 11  - 35  = 53 

min  (U,U’3  = 35 

C,  /o  ^ = 19 

(a/2 ,n^ 

^ (1 -a/2  ,n^  ~ 


Hence  C,  ) < min  (U,U’)  < C 

(a/2,nj^,n2^  ^ 

and  H^  is  accepted. 


(l-a/2 ,n^ .n^} 


(l-a/2 ,n^ ,n^) 
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12:00 


MONDAY  THURSDAY 

10:30  am  10:30  am 


No.  o£  Passenger  loading  No.  of  Passenger  Loading 

n^=8  n2=ll 

• Trains  other  than  Riverside  Trains 


FIGURE  6.  DISTRIBUTION  OF  LOADING  PASSENGERS:  CASE  3 
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MONDAY 
10:30  am 


THURSDAY 
10:30  am 


No.  o£  Passenger  Unloading  No,  of  Passenger  Unloading 


• Trains  other  than  Riverside  Trains 


FIGURE  7. 


DISTRIBUTION  OF  UNLOADING  PASSENGERS: 


CASE  3 


V.  Statistical  Model  to  Predict  Data  Gaps 
A.  Background 

Tables  4 through  7 show  four  two  dimensional  matrices  with 
"station"  as  columns  and  "time  period"  as  rows.  The  half  hourly 
periods  are  used  for  two  reasons.  First,  a uniform  time  scale  is 
required  for  all  stations  if  all  the  information  is  pooled 
together  in  the  statistical  m.odel  such  that  a prediction  of 
passenger  demand  at  Park  Street,  for  example,  will  depend  on 
the  knowledge  of  passenger  demand  for  the  same  time  span  at 
other  stations.  Secondly,  because  of  the  variable  headways 
of  the  train  arrival,  a thirty-minute  interval  is  the 
smallest  time  interval  for  which  a meaningful  estimate  of 
average  passenger  arrival  rate  can  be  obtained.  Cutting  the 
interval  into  anything  smaller  would  introduce  too  many  data 
gaps . 

Hence  Tables  4 through  7 in  fact  present  four  matrices  con- 
taining the  average  passenger  rate  of  flow  for  each  time  per- 
iod and  each  station,  on  each  of  which  a two-factor  analysis  of 
variance  can  be  performed,  a technique  that  investigates  the 
simultaneous  effects  of  two  factors,  namely  "station"  and 
"time"  on  the  response  cell,  the  passenger  loading  and  un- 
loading rate.  That  is,  if  we  designate  Rjj  be  the  average 
passenger  rate  for  station  J and  time  period  I,  Rjj  can  be 
expressed  as  a sum  of  five  components: 

^IJk  " P + + 3 + 


25 


where  y is  the  grand  mean,  a^,  the  effect  contributed  by  time 

period  I,  6-  the  effect  contributed  by  station  J,  the 

IJ 

effect  contributed  by  the  interaction  between  station  I and 
time  J,  and  the  random  error  term. 

Certain  constraints  are  placed  on  Model  (1).  They  are: 

I J I J 

2 

^IJK's  uncorrelated  with  mean  zero  and  variance  o . 

This  model  is  generally  called  a "fixed  effect"  two-way 
ANOVA  since  the  factor  levels  I and  J are  of  intrinsic  in- 
terest themselves.  If  all  the  constraints  are  satisfied, 

^IJK  then  independently  distributed  as 
N(y  + aj  + + (aB)jj,  o^) . 

B.  Modification  of  Model  (1) 

Model  Cl)  will  be  modified  as  follows: 

0 assumption  of  unimportant  interaction  effects,  or, 
for  all  I,J  pairs,  Ca6)jj  = 0 
o application  of  a generalized  linear  model  because  of 
missing  data  or  empty  cells. 

Table  application  of  a straightforward  two-way  ANOVA  to 
the  matrices  in  Tables  4 through  7 seems  simple  enough.  However, 
there  is  at  most  only  one  observation  per  cell  (i.e.  K = 1) , not 
to  mention  the  empty  cells  for  which  data  is  lacking  and  the 
prediction  of  which  is  precisely  the  objective  this  model 
set  out  to  fulfill.  Having  one  observation  per  cell 
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implies  the  infeasibility  of  estimating  the  error  term  with- 


significant. Empty  cells  make  it  inappropriate  to  use  any 
standard  ANOVA  statistical  computer  package  which  usually  re 
quires  balanced  data  (i.e.  equal  sample  size  per  cell). 
Setting  (a6)jj  = 0,  Model  (1)  then  becomes: 


a)  Regression  Approach  to  the  Analysis  of  Variance  Problem 
l^^hen  the  cells  of  the  ANOVA  matrix  have  unequal  sample 
sizes  or  empty  cells,  the  regular  ANOVA  equations  are  no 
longer  of  a simple  structure.  The  best  alternative  is  to  re 
express  the  ANOVA  model  by  its  regression  model  counterpart. 
In  formulating  the  regression  model,  we  use  indicator  vari- 
ables to  represent  the  ’’station"  effects  and  the  "time"  ef- 
fects. Specifically  we  have: 


out  assuming  the  interaction  effect  (aB)^^  be  zero  or  in- 


(l.a) 


+ £ 


IJ 


(2) 


1 


J 


T.  = 1 
1 


if  i = I and  I ^ t 


-1  for  all  i if  I = t 


0 otherwise 
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Sj  = 1 if  j = J and  J s 

= -1  for  all  j if  J = s 

= 0 otherwise 


for  all  i = 1,2,3,...  t-l,t  being  the  number  of  time 

periods 

and  j = 1,2,3,...  s-l,s  being  the  number  of  stations 
y,  are  parameters  to  be  estimated,  and  the  usual 

random  error  term;  are  uncorrelated  with  mean  E(eT.-r)  = 0 


and  variance  = o' 


IJ- 


We  note  that  the  regression  approach  (Model  (2))  gives 
rise  to  the  same  result  as  the  ANOVA  Model  (l.a).  The  in- 
dicator variables  are  coded  in  such  a way  that  their  regres- 
sion coefficients  correspond  directly  to  the  level  effects  in 
Model  (l.a).  To  illustrate,  for  a particular  I,J-  pair, 

I ^ t and  J 7^  s , 


Model 

(l.a)  gives : 

'^ij  = 

^ “l 

Model 

(2)  gives: 

" '’j 

I = t 

and  J = s , 

Model 

(l.a)  gives : 

Rts  = R 

+ 

* 6s 

+ P 

ts 

Model 

(2)  gives: 

R^s  = u 

- 2 

a . - 
1 

i 

j 

i = 

1,..  . 

t-1  and 

j = 

1,.  . . 

s - 1 . 
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If  one  equates  of  Model  (l.a)  to  that  of  Model  (2), 

the  set  of  simultaneous  equations  will  result  in  the  fol- 
lowing identities: 
y = u 

Uj  = aj  for  all  I 7^  t 

2 i = 1, . . .t-: 

i 

' for  all  J s 

= - I >=3  = ‘’s  j 

j 


Still  the  regression  approach  assumes  the  in 
is  insignificant  because  the  addition  of  int 
(TS) , for  all  possible  i,j  pairs,  makes  th 
cumbersome.  Moreover,  the  number  of  observa 
may  not  support  the  estimation  of  the  vast  n 
meters  (or  regression  coefficients).  From  t 
surface  station  data,  we  are  shown  that  each 
has  a constant  market  share  of  passenger  dem 
time  period.  For  the  underground  stations, 
of  which  are  transfer  points  such  as  Park  St 
Center,  Haymarket , and  North  Station,  additi 
the  assumption  is  necessary. 


teraction  effect 
eraction  variable, 
e equation 
tions  available 
umber  of  para- 
he  analysis  on 
surface  station 
and  during  each 
however,  many 
reet  , Government 
onal  evidence  about 
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b)  Tukey  One  Degree  of  Freedom  Test  for  Interaction 

Tukey  postulated  that,  in  an  ANOVA  model,  if  (a6) 
is  any  second  degree  polynomial  function  of  ctj  and  3j>  then 
the  interaction  term  would  be  in  the  form: 

(aS) jj  = DUjBj 

for  some  constant  D,  to  satisfy  the  constraints  set  forth 
under  Model  (1).  The  interaction  sum  of  squares,  SSAB,  will 
be : 


SSAB  = 2 I =2  2 


n2  2.2 

D OjBj 


J 


J 


Assuming  a^,  Bj  are  known,  the  least  square  estimator  of  D is 


2 2 


D = - 


I J 


2 "I  I 


I = l,...t,  J=l,...s. 


J 


The  test  statistic  for  interaction,  FSTAT,  when  only  one 
observation  per  cell  is  available,  is  then  a function  of 
the  approximate  sum  of  squares  for  interaction  (SSAB)  and  the 
sum  of  squares  of  error  (SSE) . 


FSTAT  = 


SSAB  / (SSE  - SSAB) 


1 / (df  for  SSE  - 1) 
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It  is  important  in  our  estimation  of  the  distribution  of  the 
passenger  demand  to  assure  the  correctness  of  all  the 
assumptions  we  have  set  forth  in  the  beginning.  Our  hypothe- 
sis is  : 

= Constant  market  shares  for  underground  stations, 
no  interaction  between  "station"  and  "time". 

■=  Interaction  exists,  A particular  station  at  a 
particular  time  may  have  unusually  heavy  or  light 
traffic  than  what  other  stations  normally  call  for. 
The  test  statistic  FSTAT  is  distributed  F(l,  df(SSE)-l),  and 
large  values  of  FSTAT  would  lead  to  conclusion  H^. 

C.  Presentation  of  Results 

It  is  intended  in  this  section  to  dwell  on  the  regression 
results  of  the  input  matrix  from  Table  4 in  detail 
and  briefly  discuss  the  others  since  similar  analyses  are 
performed  on  each  of  them.  Table  8 presents  the  regression 
results  on  RIVOM.MTY.  There  are  t = 25  tim.e  periods  and 
s = 11  stations,  making  the  number  of  parameters  (25-1)  + 
(11-1)  = 34  plus  an  intercept  term.  The  total  number  of 
observations  (loading  passenger  rate)  is  n=120. 

An  examination  of  Table  8 reveals  the  following  facts 
about  the  aptness'  of  the  Model  (2)  when  applied  to  this 
set  of  data: 

0 The  regression  shows  an  F value  of  6.78  which,  when 
compared  to  an  f(24.85)  at  a signficance  level  of 
a = .05,  is  highly  significant.  The  multiple 
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TABLE  8.  ANOVA  RESULTS 

ON 

RIVON.MTY: 

MODEL  (2)  ANALYSIS  OF 

VARIANCE  FOR 

THE 

REGRESSION 

(UNTRANSFORMED 

DATAJ 

SOUKCl  OF  VAlUATION 

OF 

SUM  OF  SQUARE 

MEAN  SQUARE 

F-VALUE 

ATTHIf  UTAHLE  TO  REC.I>'ESS  101^ 

^4, 

1028,2714 

30,2433 

6,7834 

DEVIATION  FKUM  kEGRESSION 

»5, 

378,965b 

4,4584 

INTEKCI-PTs  2,975  o 

multiple  CUHK.  COt;FF,»,B548l  R --  .73 
STAnL’AI-'U  EKMOrt  OF  ESTIMATES  2,111 


VApI AhLE 

REGR,  CUEFF, 

STD, error, 
OF  COEFF, 

COMPUTED 

t-value 

i 

-1  ,9401 

2,1768 

-0,8913 

3 

-2,2433 

1,2277 

•1,8273 

4 

0,7824 

0,7b92 

1,0173 

5 

1 ,6262 

0,7b92 

2,1143 

0 

-0,4358 

0,8249 

•0,5283 

7 

-0,8862 

0,8929 

•0,9925 

H 

-1,2497 

0,8169 

-1.5298 

9 

-l,l73l 

0,8204 

•1,4300 

1 0 

-0,4906 

0,9767 

-0,5105 

1 1 

-0,3600 

1,0977 

-0,3279 

1 2 

-0,3913 

1,2728 

•0,3074 

1 3 

0.2282 

1,2353 

0,1847 

1 4 

-1,7326 

1.2827 

•1,3508 

15 

-1,7481 

1,0932 

-1,5991 

lo 

-0,7497 

0,8170 

•0,9177 

17 

-1 ,3140 

0,8170 

-1,6084 

14 

-0,2261 

0,9615 

•0,2352 

19 

1.0489 

1,0724 

0,9781 

20 

0.7936 

0,9520 

0,8336 

21 

0,151t> 

0,9520 

0,1592 

22 

2.9096 

0.9520 

3.0563 

23 

4,7028 

1.0614 

4.4305 

24 

1,7016 

0,9520 

1,7874 

25 

1,4157 

1,2179 

1,1624 

2t) 

•2,3406 

0,4822 

•4,8540 

27 

-1.7365 

0,9928 

-1,7491 

28 

-0,4649 

0,7845 

-0,5926 

29 

-1,6041 

0.726O 

-2,2078 

30 

-2,1990 

0,6375 

•3,4492 

31 

4,6118 

0,5378 

8,5759 

32 

1,7968 

0,6360 

2,8252 

33 

•0,1519 

0,5641 

•0,2693 

34 

2.7625 

0,7123 

3.8784 

35 

•3,0092 

1,0478 

•2,8718 

THE  INTEHaCTION  TEST  STaT,  , FSTATI  78,1136 


ORIG.RATe  fitted  rate  residual  stand, res. 


Of  Ou 

• 1,31 

1.31 

0,62 

O.OU 

• 0 ,70 

0,7u 

0,33 

0,57 

0,57 

0,00 

0.00 

0,00 

-0,57 

0,57 

0,27 
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correlation  coefficient  = .73  indicates  that  the 
regression  has  explained  about  73%  of  the  variation  in 
the  passenger  data. 

o Small  t values  are  attached  to  some  of  the  regression 
coefficients  showing  some  stations  or  time  periods  have 
similar  mean  levels  of  passenger  traffic.  The  usual 
approach  of  eliminating  an  insignificant  variable  does 
not  apply  when  the  objective  is  not  determining  the 
correlation  between  variables,  but  rather  the  estima- 
tion of  level  means. 

o The  FSTAT  statistic  for  testing  interaction  between 
the  station  factor  and  the  time  factor  is  disturbingly 
large  (FSAT  = 78.11).  At  any  significance  level,  it 
is  concluded  that  our  assumption  of  = 0 in 

Model  (l.a)  is  not  justified.  To  remove  the  existing 
interaction,  a transformation  of  the  data  may  prove 
useful . 

By  examining  the  residuals,  additional  information 

can  be  obtained  and  the  aptness  of  the  model  determined. 

0 Figure  8 shows  a plot  of  the  residuals,  £-r-r,  vs.  the 

X «J 

/*V 

fitted  values,  Rjj-  It  reveals  a gradual  increase 
in  the  variation  of  the  residuals  around  the  expected 

/N 

mean  E(Cjj)  = 0.  As  the  value  of  R^j  increases,  so 
also  does  the  deviation  from  the  mean.  Heteroschedasity , 
or  unequal  variance  in  the  residuals,  is  apparent. 
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FIGURE  8.  PLOT  OF  THE  RESIDUALS,  E^x,  VS.  THE  PREDICTED  VALUE 


o Figure  9 shows  a normal  probability  plot  of  in 

which  the  cumulative  distribution  of  the  standard- 
ized residuals  is  compared  against  that  of  a standard 
normal  dist irbut ion , N(  0,1).  Again,  it  shows  a 
slight  departure  from  the  normality  assumption  near 
the  lower  end  of  the  scale.  As  normality  is  a rather 
robust  condition,  a slight  departure  will  not  hinder 
the  effectiveness  of  the  model. 

/*S 

Figure  10  plots  the  fitted  values,  Rjj»  against  the 
original  values  ^^jj-  This  is  done  in  an  attempt  to  detect  any 
systematic  departure  of  the  existing  model  from  the  true  model. 

/N 

A random  distribution  of  the  points  around  the  line  Rjj  = R^j 
will  give  credence  that  the  regression  model  (2)  is  indeed 
suitable.  We  notice,  however,  that  for  the  lower  values,  the 

/*v 

R 's  are  usually  greater  than  the  Rj-r's  and  the  pattern 

1 J i 

is  reverse  for  the  higher  values.  We  are  led  to  believe  then 

/N 

that  Rjj  is  some  function  of  R^j  and  a transformation  of  the 
data  is  necessary  to  straighten  the  curve. 

Transformation  of  data  to  improve  Model  (2) 

Based  on  the  information  gleaned  from  the  residual  plots 
and  ANOVA  results  above,  a modification  of  the  Model  (2)  is 
necessary.  Transformation  of  the  data  serves  several 
ob j ect ives : 

o to  stabilize  the  variance  of  the  residuals, 

0 to  improve  the  normality  of  the  residuals, 
o to  improve  the  fitness  of  tlie  model, 
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FIGURE  9.  NORMALITY  PLOT  OF  THE  STANDAPDIZED 
RESIDUALS  (UNTRANSFORMED  DATA) 
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FIGURE  10.  PLOT  OF  PREDICTED  VALUES,  R, , , VS.  THE  ACTUAL 
(UNTRANSFORMED  DATA) 


o 


to  remove  or  minimize  the  interaction  effect 


between  the  station  factor  and  the  time 
factor . 

By  examining  Figure  10,  a logarithmic  transformation 
Rjj  = log^Rjj,  for  all  pairs  seems  appropriate.  Model  (2) 
thus  becomes: 

log  R = u + y a.T.  + / b.S.  + e,  , (3) 

*e  IJ  1 1 1 J Id 

i j 


Or,  equivalently 


R 


Id 


T S 

^ n V c.  i d.  ^ ^ 

11  1 1 


Id 


where  u 

^Id 

i 

I 


log  V,  a^  = log  c^,  b^  = log  d^ 
log  njj 

1.0.  .t-1;  j = l,2,...,s-l 

1.. ..t  J = l,...s 


It  is  gratifying  to  see  the  improvements  in  the  subse- 
quent regression  results  of  Model  (3),  which  are  presented 
in  Table  9.  The  F statistic  (which  indicates  the  statistical 
significance  of  the  model)  has  increased  from  6.78  to  11.83. 
The  multiple  correlation  coefficient  square  shows  that  the 
new  regression  explains  about  83%  (compared  to  73%)  of  the 
data  variation.  The  FSTAT  has  dramatically  fallen  from 
78.11  to  a low  4.09.  Moreover,  the  residual  plots  (Figures 
11  - 13)  all  verify  the  assumption  that  Model  (3)  is  based 
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TABLE  9. 

ANOVA  RESULTS 

OF 

RIVON.MTY: 

MODEL  C3)  ANALYSIS  OF 

VARIANCE  FOR 

THE 

REGRESSION 

(TRANSFORMED 

DATA) 

SOuRCk  OF 

VARIATION 

OF 

SUN  OF  square 

mean  SOUARF 

F-VALUE 

ATTRIrtUTABLE  TO  HEGRESSION 

34, 

123,3202 

3,6273 

11,8273 

deviation 

From  regression 

05, 

26,0006 

0,3067 

INfERCI-PTa  C.481 

MULTlPl'E  CORK.  COEFF,  SOU  ABEs  , 8 255  I 
STANPAFP  EPPOR  OF  ESTIMAIE*  0,554 


VArI AhLE 

REGR,  COEFF, 

STD, ERROR, 

COMPUTED 

or  COEFF, 

T-VALUE 

2 

-1 ,1209 

0,5709 

-1,9634 

3 

-0,5345 

0,3220 

-1,6599 

4 

0,2494 

0.2017 

1,2364 

5 

0,1868 

0,2017 

0,9260 

6 

-0,0201 

0,2163 

-0,0929 

7 

-0,4322 

0,2342 

-1.8456 

R 

-0,1061 

0,2143 

•0.8686 

9 

-0 ,3063 

0,2152 

•1,7954 

lU 

• f> , 6 8 7 0 

0,2562 

•2,6940 

1 1 

-0,3178 

0,2879 

-1,1040 

1 2 

-0,2577 

0,3336 

-0,7719 

13 

0,1790 

0,3240 

0,5550 

1 4 

-0 , 2o97 

0,3364 

-0,6234 

15 

-0,3449 

0,2867 

-1,2026 

10 

-0,0106 

0,2143 

•0,0066 

17 

-0,1063 

0,2143 

•0,9696 

IR 

0,1200 

0,2522 

0,5076 

19 

0, 1 82o 

0,2813 

0,6471 

2>' 

0,4810 

0,2497 

1,9262 

21 

0,2191 

0,2497 

0,8775 

22 

0,7460 

0,2497 

2,9070 

23 

1,0444 

0,2794 

3,7517 

24 

0,6616 

0,2497 

2,6497 

25 

0,4111 

0,3194 

1,2869 

2b 

-0,0815 

0,1265 

•6,9701 

27 

-0,9117 

0,2bO4 

•3,5014 

2H 

0,0700 

0,2057 

0,3791 

2'i 

-0 , 357o 

0,1906 

-1,8732 

3C 

-0,5403 

0,1672 

-3.2313 

31 

1 ,3248 

0,1410 

9,3927 

32 

0,7701 

0,lb68 

4,6645 

33 

0,4194 

0.1479 

2,8344 

34 

0,9577 

0,1868 

5,1263 

35 

-1,9443 

0,2748 

-7,0749 

THE  interaction  TEST 

STaT,  , 

FSTATI  4,0902 

ORIG. 

FITTED 

RESIDUAL  STAND. 

RATE 

RATE 

RES. 

0 , C 0 

0,22 

0,00 

0.00 

0,00 

0,21 

0,00 

0,00 

0.57 

0,57 

0,00 

0,00 

O.O-' 

0.37 

0,00 

0,00 

39 


TABLE  9.  (CONT.) 


.RATE 

FITTED  RATE 

RESIDUAL 

STAND. RES 

oc 

7,20 

0,00 

0,00 

00 

4.22 

0,00 

0,00 

00 

2.94 

0,00 

0,00 

00 

5,04 

0,00 

0,00 

00 

0,28 

0,00 

0,00 

00 

5,68 

0,00 

0,00 

00 

0,54 

0,00 

0,00 

00 

0,53 

0,00 

0,00 

71 

i.42 

0,19 

0.34 

00 

0,92 

0,00 

0,00 

00 

0,76 

0,00 

0.00 

00 

4.93 

0,00 

0,00 

00 

2,86 

0,00 

0,00 

00 

1,99 

0,00 

0,00 

94 

3,42 

-0,15 

-0,27 

00 

0,19 

0,00 

0,00 

71 

3,85 

-0,04 

•0 ,07 

05 

0,47 

0,79 

1.43 

00 

0.46 

0,00 

0,00 

57 

1.24 

0,24 

0,43 

00 

0.80 

0,00 

0,00 

00 

0.67 

0,00 

0,00 

00 

4.31 

0,00 

0,00 

00 

2.49 

0,00 

0,00 

00 

1.74 

0,00 

0,00 

00 

2.99 

-1,09 

•1,97 

00 

0.16 

0,00 

0,00 

5B 

3.36 

0,06 

0.11 

58 

0.66 

0,88 

1,58 

00 

0.64 

0,00 

0,00 

40 

1.72 

0,70 

1.27 

00 

l.ll 

0,00 

0,00 

94 

0,92 

0,U2 

0,03 

24 

5.97 

-0,13 

•0,24 

00 

3.46 

0,00 

0,00 

00 

2.41 

-0,88 

-1.59 

75 

4,14 

-0.86 

-1,55 

00 

0,23 

0 , Ou 

0,00 

Jb 

4,66 

0,28 

0,50 

69 

0,56 

0.22 

0,39 

00 

0,54 

0,00 

0,00 

00 

1,45 

0,32 

0,58 

O'w 

0,94 

0,00 

0,00 

77 

0,78 

-0,02 

•0,03 

55 

5,05 

-0,10 

•0,19 

00 

2.92 

0,00 

0,00 

93 

2,04 

-0,06 

•0,10 

77 

3,50 

-0.23 

•0,42 

0 0 

0,19 

0,00 

0,00 

47 

3,94 

-0,13 

-0,23 

92 

0,76 

0,19 

0,34 

00 

0,74 

0,00 

0,00 

00 

1,99 

0,00 

0,00 

00 

1,29 

0,00 

0,00 

21 

1 .07 

0,12 

0,22 

9o 

6,91 

-0,00 

-0,00 

00 

4,00 

0 , Ou 

0,00 

45 

2,«o 

-0,13 

-0.24 

0 0 

4,79 

0,00 

0,00 

O') 

0,26 

0,00 

0,00 

52 

5,40 

-0,18 

-0.32 

40 
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ORIG.RATE 

FITTED  RATE 

OtOC 

0,31 

O,0‘; 

1,98 

0.0'j 

1,15 

0.00 

0,80 

0.00 

1,37 

0.0'J 

0,00 

0.00 

1 ,55 

0,00 

0,39 

n.O^ 

0,38 

0,73 

1 ,02 

O.Ov 

0,66 

0,00 

0,55 

3,22 

3,56 

o.oo 

2,06 

2.24 

1,44 

0,0  0 

2,47 

0.00 

0,14 

0,00 

2,78 

I .0  0 

0,86 

0.00 

0,83 

0.8i 

2,24 

3.11 

1,45 

0.00 

1,21 

9.53 

7,81 

1.35 

4,52 

4.87 

3,16 

11.68 

5,41 

0.30 

0,30 

0.00 

6,09 

0.19 

0,81 

0.00 

0,78 

1.14 

2,11 

2.39 

1,36 

0.00 

1,14 

13.41 

7,33 

5.75 

4,24 

2.06 

2.97 

13.19 

5,08 

0.28 

0,28 

0.00 

5,72 

0.53 

0,66 

0.00 

0,64 

0.00 

1,71 

0.84 

1,11 

0.00 

0,92 

6.33 

5,96 

2.81 

3,45 

5.22 

2,41 

3.90 

4,13 

0,21 

0,23 

0,00 

4,65 

0,10 

0,43 

0.00 

0,42 

0,00 

1,13 

0,78 

0,73 

0,00 

0,61 

4,18 

3,95 

4,73 

2,29 

0,00 

1,60 

4,00 

2,74 

0,16 

0,15 

(CONT.) 


RESIDUAL 

STAND. RES 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0 ,0o 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

-0,34 

•0,61 

0,00 

0,00 

0,00 

0,00 

•0,10 

-0,18 

0,00 

0,00 

0,44 

0,80 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,15 

0,27 

0,00 

0,00 

-1,01 

-1,82 

0,76 

1.37 

0,00 

0,00 

0,09 

0.16 

-1,21 

-2,18 

0,43 

0,78 

0,77 

1.39 

0.01 

0.02 

0,00 

0,00 

-1,45 

-2,61 

0,00 

0,00 

•0,61 

•1,11 

0,56 

1 .01 

0,00 

0,00 

0,60 

1,09 

0.30 

0,55 

•0,36 

•0,66 

0,95 

1.72 

0,00 

0.01 

0,00 

0.00 

-0,21 

•0,39 

0,00 

0,00 

0,00 

0,00 

-0,28 

•0,50 

0,00 

0,00 

0,06 

0,11 

-0,21 

•0,37 

0,77 

1,39 

-0,06 

-0,10 

-0,08 

-0,14 

0,00 

0,00 

-1  ,•47 

-2,65 

0,00 

0,00 

0,00 

0,00 

0,06 

0,11 

0,00 

0,00 

0,06 

0,10 

0,73 

1.31 

0,00 

0,00 

0,56 

1,02 

0,06 

0,11 

41 
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ORIG.RATE 

FITTED  RATE 

RESIDUAL 

STAND. RES. 

0*B3 

0,80 

0,03 

0,06 

OtO" 

O.^fB 

0,00 

0,06 

o.oy 

2,10 

0,00 

0,00 

o.oc 

1,36 

0,00 

0,00 

0 .O'"! 

1.13 

0 , Ov 

0.00 

17.33 

7,30 

0,52 

0,95 

0.00 

4,22 

0,00 

0,00 

1 .69 

2,95 

-0,5b 

•1 , 0 1 

0.00 

5,0b 

0,00 

0,00 

0.00 

0,28 

0,00 

0,00 

5.70 

5,70 

0,00 

0,00 

1 .21 

1 ,08 

0,11 

0,20 

O.Oo 

1 ,05 

0 • Ou 

0,00 

0.00 

2,83 

0,00 

0,00 

0.00 

1 ,«3 

0 ,ou 

0,00 

1.80 

1 ,52 

0,17 

0,30 

5.H 

9,84 

-0,65 

-1 .17 

8.57 

5,70 

0,41 

0,74 

3.84 

3,98 

-0,04 

•0,06 

0.0  0 

6,82 

0,00 

0,00 

0.00 

0,37 

0,00 

0,00 

o.ou 

7,68 

0,00 

0,00 

1 .00 

0,83 

0,24 

0,43 

0.00 

0,81 

0,00 

0,00 

O.Oo 

2,18 

0,00 

0,00 

0.00 

1.41 

0,00 

0,00 

1.28 

1,17 

0,09 

0,16 

7.11 

7,57 

-0,06 

•0,11 

6.25 

4,38 

0,35 

0,64 

1.65 

3,06 

-0,62 

-1.12 

0.00 

5,25 

0,00 

0,00 

0.00 

0,29 

0,00 

0,00 

0.00 

5.91 

0,00 

0,00 

2.15 

1.41 

0,42 

0,76 

0.00 

1.37 

0,00 

0,00 

o.oo 

3,69 

0,00 

0,00 

O.OO 

2,39 

0,00 

0,00 

2.10 

1,99 

0,06 

0,10 

15.56 

12.83 

0,19 

0,35 

9.18 

7.43 

0.21 

0,38 

2.15 

5,19 

• 0,88 

•1,59 

0.00 

8,89 

0,00 

0,00 

0.00 

0,49 

0,00 

0,00 

0.00 

10,01 

0,00 

0,00 

2.9b 

1,90 

0,44 

0,80 

0,00 

1,85 

0,00 

0,00 

0.00 

4,97 

0,00 

0,00 

0.00 

3.22 

0,00 

o,00 

1.52 

2.68 

-0,57 

•1,02 

15.23 

17,29 

•0,13 

-0,23 

12.87 

10,01 

0,25 

0,45 

0.00 

6,99 

0,00 

0,00 

0.00 

11,98 

0,00 

0,00 

0.00 

0.66 

O.Oo 

0,00 

0.00 

13.49 

0,00 

0,00 

1.65 

1,30 

0,24 

0,43 

0.00 

1,26 

0,00 

0,00 

0.00 

3.39 

0,00 

0,00 

0.00 

2,19 

0,00 

0,00 

1 .50 

1,83 

-0,2U 

-0 , 36 

9.23 

11.79 

-0,24 

•0,44 

42 
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• RATE 

FITTED  RATE 

RESIDUAL 

STAND. RES 

00 

3.08 

0,00 

0,00 

57 

0,5b 

0,02 

0,04 

00 

0,54 

0,62 

1.11 

00 

1,45 

0,00 

0,00 

00 

0,94 

0,00 

0,00 

54 

0,78 

-0,37 

-0,b7 

50 

5,05 

-0,37 

-0,66 

56 

2,92 

-0,62 

•1.11 

73 

2,04 

0,6o 

1,09 

90 

3,50 

0,11 

0,20 

00 

0,19 

0,00 

0,00 

00 

3,94 

0,00 

0,00 

36 

0,46 

-0,23 

•0,42 

38 

0,44 

-0,15 

•0,27 

CO 

1.19 

0,00 

0,00 

15 

0,77 

0,40 

0,73 

60 

0,64 

-0,07 

•0.12 

57 

4,13 

0,10 

0,l8 

71 

2,39 

-0,34 

•0,61 

22 

1,67 

0,28 

0,51 

00 

2.8b 

0,00 

0,00 

00 

0,16 

0,00 

0,00 

0 0 

3,23 

0,00 

0,00 

39 

0,34 

0,15 

0,27 

2b 

0,33 

-0,23 

•0,41 

00 

0,88 

0,0o 

0,00 

20 

0.57 

-1  ,05 

•1,89 

63 

0,47 

0,29 

0,51 

OG 

3,0b 

0,00 

0,00 

0" 

1.77 

0,00 

0,00 

87 

1.24 

0,84 

1.52 

00 

2,12 

0,00 

0,00 

00 

0.12 

0,00 

0,00 

00 

2,39 

0,00 

0,00 

48 

0,49 

-0,02 

•0 » 03 

52 

0,47 

0,09 

0,17 

00 

1.27 

0,00 

0,00 

47 

0,82 

-0,5b 

•1  »oi 

1 1 

0,69 

0,40 

0,87 

Ow 

4,43 

0,00 

0,00 

00 

2.56 

0,00 

0,00 

00 

1.79 

0,00 

0,00 

ui 

3,07 

0,0O 

0,00 

Of 

0,17 

0,00 

0,00 

00 

3,4b 

0,00 

0,00 

5 4 

0,52 

0,04 

0 , u9 

3b 

0,50 

-0,33 

•0  , bO 

OG 

1.35 

0,00 

0,00 

17 

0.87 

0,29 

u,53 

0 0 

0,73 

0,00 

0,00 

O': 

4,7u 

0,0o 

0,u0 

00 

2.72 

'7 ,00 

0,00 

oo 

1,90 

0,00 

0,00 

OG 

3,26 

0,00 

0,00 

O') 

0,18 

0,00 

0,00 

O'v 

3,67 

0 , Ou 

O.O'I 

5b 

0 , 9o 

— 0,32 

-0.5  8 

00 

0,78 

0,00 

0.00 

5C 

2,09 

0.51 

0.93 

1 2 

1.35 

-n,  1 9 

-0.34 

CG 

1.13 

O.f'G 

0 , OU 

4 3 
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ORIG.RATE  FITTED  RATE  RESIDUAL  STAND. RES. 


7,17 

CM 

at 

• 

0,05 

0,09 

5.55 

4.77 

0,15 

0,27 

0,00 

8,17 

0,00 

0,00 

0,00 

0,45 

0,00  ■ 

0,00 

0,00 

9,20 

0,00 

0,00 

0.81 

1,01 

-0,22 

•0,40 

0,00 

0,98 

0,00 

0,00 

0,00 

2,64 

0,00 

0,00 

0,00 

1.71 

0,00 

0,00 

0.00 

1.42 

0,00 

0,00 

10, 8u 

9,18 

0,16 

0,29 

5.63 

5,31 

0,06 

0.11 

0,00 

3,71 

0,00 

0,00 

0,00 

6,3b 

0,00 

0,00 

0,00 

0,35 

0,00 

0,00 

0.00 

7,16 

0,00 

0,00 

0,00 

0,83 

0,00 

0,00 

0,00 

0,80 

0,00 

0,00 

0.00 

2,16 

0,00 

0,00 

0.00 

1 ,4U 

0,00 

0,00 

0.00 

1,17 

0,00 

0,00 

0,00 

7,53 

0,00 

0,00 

4,36 

4,36 

0,00 

0,00 

0,00 

3,05 

0,00 

0,00 

0,00 

5,22 

0,00 

0,00 

0,00 

0,29 

0,00 

0,00 

0.00 

5.88 

0,00 

0,00 

44 


45 


FIGURE  11.  PLOT  OF  THE  RESIDUALS,  E.,,  VS.  THE  PREDICTED  VALUES, 
(TRANSFORMED  DATA) 


cum  dist.  of  the 
Stand.  N(0,1) 


FIGURE  12.  NORMALITY  PLOT  OF  THE  STANDARDIZED  RESIDUALS,  Ejj 
(TRANSFORMED  DATA) 
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FIGURE  13.  PLOT  OF  THE  PREDICTED  VALUES,  R^^,  VS.  THE  ACTUAL,  R^^ 


47 


on,  assumptions  such  as  constant  variance  and  normality  of  the 
residuals.  Note  also  that  the  estimated  values  are  positive 
throughout  while  Model  (2)  produces  negative  estimates 
occasionally.  At  this  point,  Model  (3)  is  accepted  as  a 
reasonably  good  model  for  our  purposes,  and  is  applied  to 
the  other  three  matrices  in  Tables  5 through  7,  The  results 
are  presented  in  Appendix  13.  Tables  10  through  13  present  the 
final  matrices  with  all  the  empty  cells  filled  with  the  esti- 
mated values  derived  using  Model  (3). 

VI.  Evaluation  of  Survey  Results 

Model  (3)  gives  an  estimated  time  series,  R^j  of  passenger 
rate  of  flow  at  the  underground  stations.  It  postulates  that 
the  expected  value  of  each  of  these  estimated  volumes  is  a multi- 
plicative function  of  "station”  and  "time"  effects.  Hence  the 
variance  of  these  estimates  is  again  a function  of  the 
variances  and  covariances  of  these  estimated  effects  a^ , or 
b j , etc.  Since  the  regression  is  performed  on  the  transformed 
data,  it  is  necessary  to  re-express  variances  and  expected 
values  to  their  original  forms  for  more  meaningful  interpreta- 
tion. Essentially,  Model  (3)  has: 

/S  A.  /s 

log  Rjj  = u - aj  + bj 

var  (log  Rjj)  = V(u)  + V(Sj)  + V(bj)  + 2 Cov(u,3j) 

A ^ A 

+ 2 Cov(aj,bj)  + 2 Cov(u,bj) 
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TIME 

700 

730 

800 

830 

900 

930 

1000 

1030 

1100 

1130 

J200 

1230 

l30U 

|330 

1400 

1430 

1500 

1530 

l600 

1630 

1700 

1730 

1800 

1830 

l900 


TABLE  10.  ESTIMATED  LOADING  RATE  (ALL 
RIVERSIDE  TRAINS) 


KEN. 

AUD. 

COP. 

ARL. 

BOYL. 

PARK 

GOV. 

HAY. 

N. 

SCI. 

LtCn. 

CTR. 

STA. 

PK. 

0.22 

0,21 

0,57 

0.37 

0,31 

1.98 

1.15 

0,8o 

1,37 

0,08 

' 1.55 

0,39 

0,38 

1,02 

0.66 

0,55 

3,56 

2.06 

1,44 

2,47 

0,14 

2,78 

0,86 

0.83 

2,24 

1.45 

1.21 

7.81 

4.52 

3,16 

5.41 

0,30 

6,09 

o.fli 

0,78 

2.11 

1,36 

1.14 

‘1.33 

4.24 

2,97 

5,08 

0,28 

5,72 

0,66 

0,64 

1.71 

1.11 

0,92 

5,96 

3,45 

2,41 

4.13 

0.23 

4,65 

0.43 

0.42 

1.13 

0,73 

0,61 

3.95 

2,29 

1,60 

2.74 

0,15 

3,08 

C.56 

0.54 

1.45 

0,94 

0,78 

5,05 

2.92 

2,04 

3,50 

0,19 

3,94 

0.46 

0.44 

1,19 

0,77 

0.64 

4.13 

2.39 

1,67 

2.86 

0,16 

3.23 

0,34 

0,33 

0,88 

0,57 

0,47 

3,06 

1.77 

1.24 

2,12 

0,12 

2,39 

0,49 

0,47 

1.27 

0,82 

0,69 

4.43 

2,56 

1.79 

3,07 

0,17 

3.46 

0,52 

0,50 

1.35 

0,87 

0,73 

4,70 

2,72 

1 ,90 

3,26 

0,18 

3,67 

0.80 

0,78 

2,09 

1,35 

1.13 

7,28 

4.22 

2,94 

5,04 

0,28 

5.68 

C,54 

0,53 

1.42 

0,92 

0,76 

4.93 

2,86 

1,99 

3.42 

0,19 

3,85 

0,47 

0,46 

1.24 

0,80 

0,67 

4.31 

2,49 

1.74 

2,99 

0,16 

3,36 

0,66 

0,64 

1 ,72 

1.11 

0,92 

5.97 

3,46 

2.41 

4,14 

0,23 

4,66 

0,56 

0,54 

1,45 

0,94 

0,78 

5.05 

2,92 

2,04 

3,50 

0,19 

3,94 

0,76 

0,74 

1 ,99 

1,29 

1,07 

6,91 

4,00 

2.80 

4,79 

0,26 

5,40 

o.flo 

0,78 

2,10 

1,36 

1.13 

7.3o 

4,22 

2.95 

5,06 

0,28 

5,70 

1 ,U8 

1,05 

2.83 

1.83 

1.52 

9.84 

5,70 

3.98 

6.82 

0,37 

7,68 

0,83 

0,81 

2,18 

1.41 

1.17 

7,57 

4,38 

3.06 

5.25 

0,29 

5,91 

1.41 

1.37 

3,69 

2,39 

1,99 

12.83 

7,43 

5,19 

8.89 

0,49 

10,01 

1 ,9o 

1,85 

4.97 

3.22 

2,68 

17,29 

10,01 

6,99 

11,98 

0,66 

13,49 

l,3o 

1,26 

3,39 

2,19 

1,83 

11,79 

6.82 

4.77 

8,17 

0,45 

9,20 

1,01 

0,98 

2,64 

1.71 

1.42 

9,18 

5.31 

3.71 

6,36 

0,35 

7,16 

0,83 

0,8o 

2,16 

1 .40 

1.17 

7,53 

4,36 

3,05 

5,22 

0,29 

5,88 

TABLE 

11. 

ESTIMATED  UNLOADING 

RATES 

^ (ALL 

- 

RIVERSIDE  TRAINS) 

TIME 

KEN. 

AUD. 

COP. 

ARL. 

BOYL. 

PARK 

GOV. 

HAY. 

CTR. 

700 

2,48 

1,20 

4,57 

2,82 

0,50 

2,49 

0.52 

0,08 

730 

1,38 

0,66 

3.1b 

1 ,56 

0,28 

1.22 

C,29 

0,04 

800 

3,00 

1.72 

6,68 

8,20 

0,71 

1.97 

1,00 

0.09 

830 

4,19 

2,27 

8,19 

7,91 

0.94 

4,74 

1 ,96 

0,06 

900 

5,05 

1.77 

6.77 

4,26 

0,74 

3.20 

0,57 

0,13 

930 

1,79 

0.91 

3.47 

2.17 

0,38 

1,76 

0,44 

0 , 06 

1000 

1.33 

4.33 

6,64 

4.09 

0,81 

3.25 

0,81 

0,12 

l03o 

1,44 

0,6e 

3.87 

2,65 

1,65 

l.'O 

C , 43 

0,07 

1 100 

1,78 

0,35 

1,64 

0,55 

C,2o 

0.89 

0,19 

0,03 

1 1 30 

1.55 

C,7l 

3.01 

0,80 

0,89 

1,64 

0,34 

0 ,05 

1 200 

1,68 

0,73 

3.33 

2.63 

0,36 

1.81 

0,38 

0,06 

1 230 

0,48 

0,59 

6,00 

1.35 

0,25 

1.23 

0,26 

0,04 

1 300 

1,64 

0,79 

3.03 

1,87 

0,33 

1,65 

0,34 

0,05 

1330 

1.77 

0,79 

2,80 

1,86 

0,33 

1.65 

0,34 

0,05 

1 400 

2.32 

1.43 

3.15 

3,36 

0.44 

3,48 

0,62 

0,24 

1430 

6.42 

0,98 

2.20 

2,30 

0,18 

2,28 

0.42 

0 . 0 7 

l500 

3.81 

0,88 

3,38 

2,08 

0,34 

1,87 

0,38 

0,03 

l53o 

1.07 

1,70 

6,50 

4,01 

0,20 

7,00 

0,74 

0,38 

1600 

2.6b 

2,u4 

7,78 

4,8o 

2,60 

3.38 

0,57 

0,14 

l630 

3,38 

2,47 

9,44 

5,82 

1,09 

4.31 

1,82 

0.16 

l700 

7,38 

2.34 

8,95 

5,52 

1.34 

8 , 7o 

1 ,05 

0,04 

1730 

10,04 

1.02 

6.94 

4,28 

0,19 

5.37 

0,87 

0,12 

1 800 

8,06 

2,84 

10.86 

6,69 

0,94 

4,00 

0,63 

0,5o 

l830 

2,75 

1.55 

5.91 

3,64 

0,64 

5,70 

0,44 

0,10 

1900 

0,34 

0,16 

0.62 

0,38 

0,07 

0,34 

0,07 

0,01 
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TABLE  12.  ESTIMATED  LOADING  RATES 

CALL  NORTH  STATION  TRAINS) 


TIME 

KEN. 

AUD. 

COP. 

ARL. 

BOYL. 

PARK 

GOV. 

CTR. 

HAY. 

700 

0.32 

0.23 

0.33 

0,19 

0,04 

0.23 

0,05 

0,02 

730 

0,45 

0,33 

0,47 

0,27 

0,06 

0,33 

0,07 

0.03 

800 

2.63 

1.94 

1.25 

1,56 

0,36 

8.10 

0,14 

0,17 

830 

3,06 

1,80 

3,3b 

1,84 

0,33 

2,28 

0,13 

0,16 

900 

5,81 

2,01 

2.82 

2.17 

0,37 

0,63 

0,43 

0,18 

930 

J.J4 

1.08 

1,51 

0,42 

0,20 

0,79 

0,54 

0.13 

1000 

0,35 

0,96 

1,18 

0,68 

0,33 

1.21 

0,15 

0,07 

1030 

1.14 

1,66 

1,34 

0,92 

0,07 

1,00 

0.33 

0.06 

1 100 

3,29 

1.13 

0,94 

1,00 

0.16 

0.81 

0,94 

0,12 

1130 

0,93 

1,45 

1,26 

0,36 

0,62 

0,95 

0,11 

0,08 

1200 

1,75 

0.62 

0,46 

0,56 

0,11 

0.62 

0,10 

0,05 

1230 

1,66 

0,89 

5,23 

0,71 

0,16 

0,88 

0,03 

0,08 

1300 

1.11 

0.82 

1,15 

0,66 

0,15 

0,81 

0,16 

0,07 

1330 

2.94 

1,25 

1,60 

1,01 

0,23 

1.24 

0,24 

0,07 

1400 

3. .8  4 

1,30 

2,38 

1,04 

0,24 

1,29 

0,25 

0,04 

1430 

3,9b 

1.27 

1,74 

1,02 

0.23 

1,26 

0,25 

0,05 

1500 

3,00 

2,21 

3.11 

1,78 

0.41 

2,19 

0.43 

0,19 

1530 

6.84 

5,05 

7,10 

4.05 

0,93 

5,00 

0,98 

0,44 

l600 

2.00 

2.62 

3.68 

2,10 

0.48 

1 ,8o 

1,30 

0,23 

1630 

3,67 

2,54 

4,54 

2,59 

0.63 

5.81 

0,24 

0,58 

1700 

2.14 

1,36 

4.93 

1,58 

1,09 

2,18 

0,90 

2,22 

1730 

4.35 

4.90 

5,67 

3.90 

0,53 

5,7o 

0,97 

0,24 

1 800 

1 .08 

1.42 

1.88 

2,50 

0,08 

1,00 

0,72 

0.12 

1830 

0,83 

1.13 

1.59 

1.19 

0,25 

1,11 

0,26 

0,10 

TABLE  13.  ESTIMATED  UNLOADING  RATES 
(ALL  NORTH  STATION  TRAINS) 


TIME 

KEN. 

AUD. 

COP. 

ARL. 

BOYL. 

PARK 

GOV. 

HAY. 

N. 

CTR. 

0,84 

STA. 

700 

0,19 

0,47 

0,47 

0,57 

0,25 

3,30 

2,04 

0,67 

73u 

0,39 

0,96 

0,95 

1.15 

0,51 

6,68 

4.12 

1,69 

1.3b 

80o 

1.42 

3,49 

2.69 

4,18 

1,85 

22,07 

20,43 

6,16 

5,15 

8 30 

0,50 

1.93 

3,96 

3.90 

1,02 

10,69 

11.65 

3.40 

1,09 

900 

0.39 

2.32 

2,30 

4,67 

1,23 

13,07 

16,97 

4,09 

3,53 

930 

0,59 

1.00 

0.99 

2,96 

0,53 

5,5b 

7.15 

0.47 

1,09 

1 OOo 

0,35 

2,67 

1 .25 

1.51 

0,89 

8,66 

6,88 

0,94 

1 ,90 

I03o 

0.47 

1.45 

1 ,02 

1 ,46 

0,52 

5,11 

2,93 

2,22 

1.47 

1100 

0,43 

1.45 

1,44 

1.57 

0,61 

7.19 

3.97 

0,49 

1 , 4u 

1 13o 

0,29 

0.68 

0,85 

0,62 

1.12 

4.88 

1.14 

1.24 

1,00 

1 200 

0,32 

0,97 

1.19 

0,96 

0,52 

6,78 

5,10 

1.72 

1,38 

l230 

0,88 

1.42 

1.15 

1,70 

0,75 

9 ,9n 

4.92 

2.51 

2,02 

1300 

C,33 

1.43 

1.41 

l,7l 

0,76 

9,92 

6.12 

2.51 

3,5b 

1330 

0,5y 

0.78 

0,68 

0.93 

0,41 

5,39 

3,33 

0,5o 

2.14 

1400 

1 ,09 

1.75 

2 , Ou 

2,10 

0,93 

12,19 

7,52 

2.39 

1,82 

l43o 

1 , 04 

1.41 

0,7  2 

1,69 

0,75 

9,82 

6,06 

2,2b 

2. 3b 

1500 

0,55 

1.35 

0 , 8o 

1 .61 

0,71 

9,38 

5,79 

3,69 

1,92 

1530 

0,6o 

1 ,48 

1.47 

1 .78 

0,79 

10.33 

6,37 

2,62 

2,11 

IbOO 

0,29 

1 ,02 

1,01 

1.23 

0,54 

8,17 

5,52 

1,81 

1,46 

1 6 3 0 

0,39 

0,85 

1.21 

1.47 

0,38 

16,22 

4,79 

3.9b 

1.74 

1700 

0,23 

0,71 

1,09 

0,5  0 

0,74 

10, 3o 

4,67 

9,04 

1,56 

1730 

0,39 

0,9o 

1 ,28 

0,94 

0,41 

12,33 

5,34 

9,05 

1,83 

l800 

o,50 

0,77 

1,00 

0,88 

0,46 

7,7o 

3,62 

3.41 

1 ,44 

1830 

1.37 

0,8  0 

0,79 

0,85 

0,25 

3,6b 

2,39 

1 .41 

1,14 

where  Var(log  Rjj)  represents  the  variance  o£  the  present 

estimate  log  •'ij  from  the  true  mean.  Since  an  actual  observed 

rate  of  flow  varies  from  the  true  mean  value  with  variance  , 

a predicted  value  of  any  individual  observation  will  have 
2 

variance  (a  + Var(log  Rjj))-  To  re-express  the  transformed 
data  to  their  original  form, 

R = 

IJ 


Var(Rjj)  = Var(log  R^j) 

and  the  variance  of  a predicted  individual  observation  is 

✓N  2 2 ^ 

(Rjj)  (a  + Var  (log  Rjj))-  Such  is  also  the  variance  of  the 
distribution  of  passenger  flow  to  be  calibrated  into  the 
simulation  model  to  generate  new  time  series.  We  choose  Park 
Street  Station  to  illustrate.  Figure  14  represents  three  time 
series,  namely  the  actual  collected  data,  the  estimated  data 
derived  from  Model  (3),  and  the  95%  confidence  bound  not  on 
the  estimated  series  but  rather  on  a future  predicted  observa- 
tion. That  is,  model  (3)  gives  only  an  estimated  average  rate 
of  flow,  when  given  the  "station"  and  "time"  are  fixed,  but  in  95 
cases  out  of  100,  an  individual  observation  will  fall  within  the 
confidence  bound.  The  wide  confidence  band  indicates  that 
any  actual  observed  series  of  passenger  flow  data  is  expected 
to  behave  rather  sporadically  around  the  means. 

VI I . Conclusion 

In  improving  the  operational  performance  of  the  Light 
Rail  Vehicles  in  terms  of  maximizing  the  utility  of  the  LRV 
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Pass . /Min , 
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FIGURE  14.  PROFILE  OF  LOADING  PASSENGERS  AT  PARK  STREET  STATION,  1978 
(TO  RIVERSIDE) 


cars  and  reducing  passengers’  waiting  time,  information  about 
the  present  distribution  of  ridership  of  a transit  line  is 
essential.  This  study  has  shown  that  sufficiently  accurate  and 
quantitative  knowledge  about  such  distributions  at  different 
stations  can  be  gained  without  expending  huge  amounts  of  labor  and 
financial  resources.  Through  the  use  of  sampling  techniques, 
from  sampling  days/shifts  and  stations  to  sampling  cars  or  doors 
of  cars,  and  through  the  application  of  a simple  statistical 
model  in  which  "station”  and  "time"  are  the  determining 
factors,  a representative  daily  ridership  profile  for  each 
station  has  been  obtained.  However,  there  is  still  room  for 
improvement.  For  example,  the  problem  of  sample  size  has  not 
been  addressed.  Had  a replicated  sampling  scheme  been  fol- 
lowed, that  is,  if  for  each  station-time  period  combination, 
more  than  one  observation  of  the  passenger  rate  of  flow  were 

taken,  it  would  be  possible  to  estimate  the  magnitude  of  the 

2 

interaction  effect  and  the  pure  random  error,  a , separately. 

Also,  a cursory  look  at  Table  11,  the  estimated  unloading 
passenger  rate  for  the  Riverside  Train  reveals  abnormally  high 
estimates  for  Copley  during  the  afternoon  rush  hours.  The 
environment  (mostly  business  establishments)  around  the 
locality  does  not  command  such  high  return  rates  of  passengers. 

Validation  plans  are  underway  to  verify  or  modify  our  cur- 
rent estimates.  Additional  data  for  the  12:00  - 1:30  period 
needs  to  be  gathered.  However,  instead  of  sampling  all 
stations,  it  is  proposed  that  only  three  to  four  stations 
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should  be  observed,  and  a continuous  time  series  of  passenger 
data  be  collected.  We  also  need  to  examine  which  station 
and  time  period  in  which  interaction  plays  an  important  part. 
Also,  it  is  discovered  from  the  data  sheets  that  some  sur- 
veyors often  do  not  systematically  sample  the  cars  when  count- 
ing passengers  from  all  cars  becomes  unmanageable.  This  should 
be  emphasized  to  assure  the  representativeness  of  the  col- 
lected data,  since  the  distribution  of  people  loading  in  and 
out  of  the  two  cars  is  often  not  uniform.  Depending  on  their 
proximity  to  the  subway  entrance,  one  car  usually  gets  more 
passengers  than  another. 

Finally,  a word  should  be  said  about  the  estimated  total 
volume  of  passengers.  An  earlier  study  on  surface  stations 
has  estimated  that  the  daily  total  number  of  patrons  from  the 
surface  section  of  the  Riverside  line  trains  is  only  approxi- 
mately 6,800.^  It  would  be  interesting  to  see  if  the  data  one 
year  later  (underground  stations)  in  this  study  corroborates 
previous  findings.  We  approximate  the  total  daily  volume  of 
passengers  going  outbound  to  beyond  Kenmore  as: 


Total  RIVON  outbound 


/s 


^•^RIVON 


R 


\ 


Kwok,  Betty,  "MBTA  Passenger  Demand  Analyses,  1977,” 

U.S.  Department  of  Transportation,  Urban  Mass  Trans- 
portation Administration,  UMTA-MA-06-0074-81-2 , August  1981. 
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This  is  equivalent  to  multiplying  the  ’’net  on”  rates  (pas- 
senger/min.)  by  30  (minutes) . Similarly,  the  total  number 
of  passengers  from  the  surface  stations  going  inbound  to  the 
tunnels  is  approximated  as: 


These  two  statistics  are  9,700  and  10,000  respectively.  Has 
the  ridership  really  increased  40%?  Statistical  variations 
around  earlier  estimates  and  current  estimates  prohibit  such 
a straightforward  conclusion.  We  can  only  be  contented  by 
the  apparent  correspondence  between  the  total  outbound  and 
inbound  volumes,  which  adds  to  the  credibility  and  reliability 
of  our  passenger  flow  estimates. 


TOTAL  RIVON  inbound  - 30 
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APPENDIX  1 
SURVEYOR  DATA 


DATE 


O Lechmere— > Riverside 
n^-^iverside  — > Lechmere 


Initial  of  Observer 


Remarks : 

indicate  if  count  is  not  for  entire  train 

1.  car  number  - Cl,  C2,  C3 

2.  door  number  - Dl,  D2,  D3 

indicate  any  unusual  circumstance  such  as 
big  groups,  etc. 


Time  of 
Arrival 

LINE 

Car 

Length 

PCC  / 

v/LRV 

PASS.  ON 

PASS.  OFF 

REMARKS 
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3 / 
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date  

STATION 


l/*r^echme re  — » Riverside 
□ Riverside  — ^ Lechmere 


Initial  of  Ob 


Remarks : 

indicate  if  count  is  not  for  entire  train 

1.  car  number  - Cl,  C2,  C3 

2.  door  number  - Dl,  D2,  D3 


indicate  any  unusual  circumstance  such  as 
big  groups,  etc.. 

7, 
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APPENDIX  2 
MBTA  FILE  STRUCTURE 


780522PARK 

0750H  IPCC 

780522PARK  ' 

0751BCILRV 

780522PARK 

0751H  ILRV 

TBCTSTTPA-p-K  

0752Vi  IPCC 

■ 

780b22PARK 
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20 

9 
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0 
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082^P  2L'^V 

8S 

4 

760527PAPK 

0826H  IPCC 

780522DAC  K 

0827M  1D(C 

780522PAPK 

08j1H  ILhV 

780522PARK 

0836H  l^CC 

780522PAP  K 

0837RC3PCC 

780522  PARK 

08  33CC2LP V 

780522PA3  K, 

0339CC2PCC 

780522PAOK 

0340R  2LSV 

19 

8 

730522PAPK 

0841  CCP  PCC 

730522PAP  K 

0842BC  3PCC 

760522PAC  K 

084?  BL3 

7d0522PAr  K 

06 '<3'.'  ILRV 

780522  r’A'^K 

0844H  IDCC 

7B0  522PAR  K 

0844CC2PCC 

780 52 2 PARK 

0345H  IL=V 

780522PARK 

0845CC2PCC 

780  52PPAPK 

084oH  IPfC 

7b052?PAPk 

0847P  ?LPV 

10 

i 1 

780  52?  '’APR 

08480  ?LPV 

5 

0 

7805??OA^k 

0849CC2PCC 

78052'’PAP  K 

0O51H  IPCC 

78052Pf’A''K 

08 5 IN  ILRV 

780522  PA  = R 

08  5 4k  ’LhV 

1 

1 

730522  K 

0d54H  LL"'V 

780522PA'^K 

0S5n3C2PC: 

780522  PAPK 

088^ri  IP^C 

730522PARK 

09J0H  l°CC 

780522  PARK 

OOJOH  IPCC 

780522PAP  K 

OSGl CC?  C 

780522PA^K 

0902H  l.PV 

780522  ^At:  K 

090?riC3PCC 

780522PAKK 

0O04N  ILRV 

780522  K 

0904H  locc 

780522PAP  K 

0907k  ILhV 

7805?2'5ARK 

0908CC2OCC 

7«0  522P'\k  K 

090O8C1L0V 

7805  2?  P AInK 

0O09H  3L-V 

7r)05?2?A'^ 

09  10  8 2 LPV 

13 

4 

X X 


X X 


X X 


X X 
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APPENDIX  3 
PROGRAM  EXTRAC 


TYPE  EXTRAC 
OPEN  MBTA. 

F STATION  EQ  "ARLINGTON". 

INI  1 ARLING.TON. 

SORT  BY  MAYDAY  LINE  TIME  TYPEOFCAR  TRANSIZE. 

PRINT  ON  1 MAYDAY  TIME  LINE  TRAINSIZE  TYPEOFCAR  PASSON  PASSOFF  $ 
CAR!  CAR2  CARS  DOORl  DOOR2  DOORS  COLLECT  FMT  D4  14  A8  II  AS  & 

IS  IS  F.l  A1  A1  A1  A1  A1  27  II  END. 

QUIT. 
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APPENDIX  4 

PROGRAMS  EDIT  AND  EDITl 


PROGRAM  EDIT 


INTtGER  LINE,  CT,  CARl,  CAR2,  CAR3,  Dl,  D2,  03 

--CAL-fc  

CALL  0F1LEC22, 'STRAI . ' ) 

J --  -fCltn-rt) 

1CU)2  = 0 

- - - ICt03=t) - 

1C?01=0 

IC7D2-0  - - - 

IC2Li  3 = 0 


1C3D2SO 

IC30  3SO  — - - 

REALC  2t  ,5,  ENOslOOHUATE,  IT  I ME  , 1. 1 NE  , lTi> , CT , I UN  , I UFF  , C AR  1 , 

2 CAn2, CAR3, 01 ,02,03, ICOL  - 

5 FURl'AT(lH,14,A2,Il,A3,2(I3),6{An,2X,Il  ) 

TrrcT.EQ, 'occ' iGo  Tu  2r>  ' — • 

1F(CAR1  .EU,'X',  AM).  CAK2.F(j.  *X*  . AND,  Dl.FU.'X', 

2 Al.'O  ,03.F0.'X'l  Gn  TO  70  

U (CAl-.t  .NE.'X'IGO  Tij  7 

IFUn  .Ea.'X'')  IClOlrl  - --  - • 

U (1’2.FU.'X')IC1D2=1 

IF(r3.EQ,'X'1IC103  = l 

II-  (L-l  .NE.'X'.ANO.D2.NE,'X*.ANn,03,NE.'X')Gn  TO  b 

GO  TO  7 ■ ■ ■“  ■ ■ 


b IClOIsl 

lC1U2=t 

rc  1 03= I 

7 rFfCAP2.NE. 'X'lGO  TO  9 

IK  (L  1 .Eu.'X'UC201  = l 
IF{P?.EG. 'X'1 If202=l 
n (i)3,FO.'X')Jr203  = i 

lF(l'r.'7rE.''X'.A''IU.D?.rJE.'X'.AMD.D3.NE.'X'')G0  TU  H 
Go  TU  9 
n TC2ri  = 1 

IC21  2 = 1 
IC7D3.=  1 

9 I F ( C A iM  . ;•!  I . ' X ' . A U . C A H 2 . M E . ' X ' . A N D . 0 1 , N E , ' X * . 

2 A.-D.n2.0F.'X'.Ar;n.n3.,'JF.'X*)GU  TH  10 

G'j  ru  12 

10  Tcini=i 

1C  1 i>2=  1 
IC  103=1 
1 C 2 1 > 1 = 1 

ic?r2=i 


U'2i  3=1 

12  1 r { ( JC101  + rc 102  + 1: 1 P3  + TC201 +TC2D2  + IC2D3) .GT.OIGO 

'M-  1 H:  IS,  lA  ) n IME,  l!)AVE 

11  FOO'-  AI  i 'CAT-I.OT  ALI.nCATF  CDUOT  TO  OHOP  AT*,  IX,  14, 

^ 1 X , ' I '.  * , 1 X , I 9 ) 

GO  TO  75 


25 

IF  ( C R 1 , 

•-'E. 

*X*1G"  1" 

27 

IF (L  1 .CO 

. 'X 

'■)lCini  = i 

1 )•■  ( : ) 2 . F U 

'110  02=1 

10(1  1 . r 

' i 

• 

* .Ao'i'.u:. 

G E . * X * 

IGn 

GU  TP  27 

2o 

I C 1 0 1 = 1 

TC  1 0 2= 1 

27 

IF(CA,;2. 

OF  . 

'X'lGO  TO 

2 9 

I o ( I.M  . Ew 

. 'X 

'1  l('20l  = l 

. 'X 

' 1i  ('  2 ■')  2 = 1 

1 i-  ( 1 1 . ' U-: 

.'X 

' o?  . 

OE. *X* 

)G0 

TTT  75  ■■ 
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PROGRAM  EDIT  (CONT.) 


Ci.  r.i  i>y 


24 

ir?i 

1 = 1 

1 1:  V 1 

2=  1 

29 

1 IM  c 

AA  3.!4f-,. 'X  ' ICO  T(i  31 

1 M. 

1 .tij.  'X ' 1 ir 3.0 1 = 1 

1 f ( 1 

2.Fu.  ' ' ur  30/  = i 

1 1-  ( 1.- 

1 .'-'I-  . 'X  ' t . 'X*  IGO 

G'  • 1 

o 3 1 

3 j I C 3 1 ' 1 = I 

IC3!)2=i 

3 1 Ji-  (CAhI  .Nt.,  ‘’X'.AM/.('AK2„Nt,  'X*,AMO.ITS,t.O,2.A.<D.ni  .we.  'X'. 

2 A a),D2.f-f  . *X*  )Gll  TO  3 3 
(•(  ,U  31 

3 3 1 C 1 1 1 = I 

1C  h-2=l 
IC?L  1=1 
IC2l'2=l 

3 1 II  (CAi  U.v.-,.  ‘’X*’  .A'' j.CAKi.l-iK.  "X®  . Am;.CAK3„''lt.  '■X'.AWU. 

2 1 rSo<:;i. . 3 , A i-l'.D  1 „ 'x  " . ANP  .[>2  'X  ^ )GU  fO  37 

Gi!  3 '.I  3 r> 

37  ICtltsl 

fC  1 i‘2=l 
I C20 I = 1 
1C2U2=1 
1 C3u I = 1 
1C  31^2=1 

3ri  I K ( I 1C  1 l;l  ■*  rc  1 P2+IC2U1  4 IC2P2  ♦•TC3D  1 •>  IC  3D2  ) ,GT.  OlGl)  m 75 

«(■  II  h,(b  , 1 4 ) n I .'4t , IDATF,  . . - 

G..  'lu  7b 

7(,  lClu3=l  ' - . 


IC2U1 =1 

75  4,i^lTE(  22,80)ruATK,ITIME,LlNt:,ITS,CT,TON,TOFF,  

2 ICini,IC1U2.TClD3,lC2Dl,IC2D2,lC7D3,IC3Dl,lC3D2,IC3D3»ICOL 

fio  F0FMAT(19,T4,A2,I1  rA3, 13,13, ‘7(in,'7X,Tl)  

GU  Til  1 

10  0 C’:0  ' -----  --  — - 
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PROGRAM  EDITl 


Program  EDITl  is  similar  to  EDIT  but  is  used 
for  trains  opening  doors  on  both  sides. 


INTeGEE  line,  CT,  CAF.1,  CAE2,  CAT.3,  D1,  D2,  D3,SIDE 
CALL  IFILE(21, 'KEMOFE.  ’ ) 

CALL  OFILE (22, • SLFAI.  ' ) 

10  IC1D1L=C 

IC1D2L=0 
IC  1D3L=0 
IC2D1L=0 
IC232L=0 
1C2D3L=0 
IC3D1L=0 
1C3D2L=0 
IC3D3L=0 
IC1D1F.  = 0 
IC1D2R=0 
IC103K=0 
IC2D1R=0 
1C2D2E=0 
IC2D3R=0 
IC3D1E=Q 
IC3D2E=C 
IC3D3E=0 

. SEA  D (2  1 ,50,  END=1  000)  I DA  IE,  IT  I. 'IE, LINE,  IIS  , CT  , ION  , 10  FF,  CAE  1 , 
2 CAR2,CAR3,D1,D2,D3,ICCL,SIDE 

50  FOF.WAT  (I3,I4,A2,I1,A3,2(I3),6(A1),2X,I1,A1) 

IF  (CT- EQ-  ’PCC)  00  TO  250 

IF  (CAEl .SQ. ' X*.  AND-  CAR2.EU.'X'.  AND.  DI.EQ.'X’. 

2 AND  .D3.E0.'X')  GO  TO  700 
IF  (CAR1.ME.  'X')  00  TO  70 
IF  (SIDE.  NE.  'L'  ) GO  TO  61 
IF  (D1 . SO. ' X' ) IC1D1I=1 
IF  (D2. £0. ' X ' ) IC1D2L=1 
IF  (D3. EQ.  'X')  IC1D3L=1 
GO  TO  70 

IF (D1 . NE.  ' X ' . AND. D2. NE.  'X' .A ND. D3. NE.  * X’ ) GO  TO  60 
GO  TO  70 

60  IC1D1L=1 
IC1D2L=1 
IC1D3L=1 
GO  TO  70 

61  IF(D1.EU. 'X')  IC1D1R=1 

IF  (D2.EQ. 'X') IClD2R=1 
IF  (D3.EQ.  'XM  IC1D3H=1 

IF  (D1. NE.  *X *. AND.D2. NE.  *X' -AND. D3. NE.  ' X' ) GO  TO  62 
GO  TO  70 

62  IC1D1K=1 
IC 1D2H=1 
IC1D3R=1 

70  IF(CAR2.NE.  *X')GO  TO  90 

IF  (SIDE.  NEN  • L' ) GO  TO  81 
IF  (D1. EO. *X')  IC2D1L=1 
IF  (D2. EO. 'X')  IC2D2L=1 
IF  (D3. SO. ' X' ) IC2D3L=1 

IF  (D1.  NE. 'X'.  AND.D2.NE.  'X’-AND.DS.  NE.  'X')  GO  TO  30 
GO  TO  90 

BO  IC2D1L=1 

IC2D2L=1 
IC2D3L=1 
GO  TO  90 

81  IF (D1.E0. 'X') IC2D1R=1 

IF (D2. £0. • X» ) IC2D2R=1 
IF(D3.EQ.  'X')  IC2D3R=1 
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PROGRAM  EDITl  (CONT.) 


IF (Dl.  NE- * X' . AND. D2. NE- VX* -A  NO. D3. NE.  • X* ) GO  TO  82 
GO  TO  90 
82  IC2D1H=1 

1C2D2R=1 
IC2D3R= 1 

90  IF (CAR1 . NE. ' X' . AND.CAE2.NE. ' X' . AND. D1. NE. • X‘ . 

2 AND.D2.NE.'X'.AND.D3.  NE. 'XM  GO  TO  100 
GO  TO  120 

100  IF  (SIDE.NE.  ’!•) GO  TO  101 
IClDlL=1 

IC1D2I.=  1 
IC1D3L=1 
1C2D1L=1 
IC2D2L=1 
IC2D3L=1 
GO  TO  120 

101  IC1D1R=1 
IC  1D2R  = 1 
IC1D3P=1 
IC2D1R=1 
IC2D2R=1 
IC2D3R=1 

120  IF(  (IC1D1L  + IC1D2L*IC1D3L+IC2D1L*IC2D2L>IC2d3L)  .GT.O.OE. 

2 (IClDlf>IClD2R<-IC1D3E«.IC2D1R*IC2D2R«’IC2D3R)  .GT.O) 

2 GO  TO  750 

KEilE  (5, 140)  ITIME.  IDATE 

140  FORMAT  ('CANNOT  ALLOCATE  COUNT  TO  DOOE  AT', IX, 14, 

2 1 X, 'ON' ,1X,I8) 

GO  TO  750 

250  IF(CAP.1oNE.  'X')  GO  TO  270 

IF  (SIDE.  NE. 'L')  GO  TO  26  5 
IF (D1 . Eg. ' X' ) IC1D1L=1 
IF  (D2. EO- 'X') IC1D2L=1 
IF  (D1- NE.  ' X ' . AND. D2. NE.  'X ' ) GO  TO  260 
GO  TO  270 
260  IC1D1L=1 

IE  1D2L=1 
GO  TO  270 

265  IF  (D1.EQ. 'X') IC1D1B=1 
IF  (D2.EU« 'X') IC1D2R=1 

IF  (D1. NE. 'X'- AND.D2. NE.  ' X' ) GO  TO  266 
GO  TO  270 

266  IC1D1R=1 
IC1D2P=1 

270  IF  (CAE2.NE.  'X') GO  TO  290 

IF  (SIDE.NE.  'L' ) GO  TO  28  1 

IF(D1,E0.  'X')  IC2D1L=1 
IF  (D2. EO.  'X')  IC2D2l=1 

IF  (D1 . NE. ' X • .A ND. D2. NE.  ' X' ) GO  TO  280 
GO  TO  290 

280  IC2D1L=1 
IC2D2L= 1 
GO  TO  290 

281  IF  (Dl. EO. 'X' ) IC2D1R=1 

IF  (D2.E0. *X’) IC2D2?=1 

IF (D1 . NE. ' X • . A ND.D2. NE. ' X' ) GO  TO  282 
GO  TO  290 

282  IC2D1R=1 
1C2D2K=1 

290  IF  (CAE3. NE.  ' X* ) GO  10  310 

IF  (SIDE.  NE.  ' L' ) GO  TO  305 
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PROGRAM  EDITl  (CONT.) 


300 

305 

306 
310 

330 

335 

340 

370 


375 


380 


700 

705 

750 

BOO 

1000 


•X' ) GO  TO  300 


•XV)GO  TO  306 


,'X*.AMD-ITS.EQ.2.AND.D1.NE.'X' 


. ' X’ . AHD. 

TO  370 


IF  (DI.SQ. 'X* J IC3D1L=1 
IF  (D2.EQ. 'X') IC3D2L=1 
IF  (D1.  ME.  *X'.  AMD.  02.  ME, 

GO  TO  310 
IC3D1L=1 
IC3D2L=1 
GO  TO.  310 

IF  (D1.  EO.  'X'  ) IC3D1R  = 1 
IF  (D2.EQ.  'X')IC3d2E=1 
IF  (D1.NE.*X'.AND.D2.NE 
GO  TO  310 
IC3D1R=1 
IC3D2R=1 

IF  (CARI.  ME.  'X'-AND.CAEI.NE. 

2 AMD-  D2.ME.  'X*)  GO  TO  330 
GO  TO  340 

IF  (SIDE- ME.  • LM  GO  TO  335 
IC1D1L=1 
IC  1D2L  = 1 
IC2D1L=1 
IC2D2L=1 
GO  TO  340 
IC  1D1R=1 
IC 1D2R=1 
IC2D1R=1 
IC2D2R=1 

IF  (CAR1-ME.  'X'.AMD-CAR2-.ME.'X'-AND-CAS3.SE. 

2 itS.E0.3.AND.D1.NE.'X'.AND.D2.NE.  'X')  GO 

GO  TO  380 

IF  (SIDE.  ME.  ' I' ) GO  TO  375 
IC1D1L=1 
IC1D2L=1 
IC2D1L=1 
IC2D2L=1 
IC3D1L=1 
IC3D2L=1 
GO  TO  380 
IC1D1E=1 
IC1D2E=1 
IC2D1?.=  1 
IC2D2R=1 
IC3D1E=1 

IF^f°(iciD1L  + IC1D2L>IC2D1L^IC2D2L+IC3DlL*IC3D2L).G' 

2 OR.  (IClDlR  + lC1D2EtlC2D1R*IC2D2E+IC3D1F>IC3D2E)  , 

2 GO  TO  750 

'.’EITS(5, 140)  ITIHE,IDATE 
GO  TO  750 

IF  (SIDE.  ME.  • L* ) GO  TO  705 

IC1D3L=1 

IC2D1L=1 

GO  TO  750 

IC1D3E=1 

k’EITE  U2.80  0)  IDATE.ITIME,  LINE,  ITS,  CT,IOM,IOFF, 

2 IC1D1L,IC1D2L,IC1D3L,IC2D1L,IC2D2L, IC2D3L,_  IC2D2E 

2 lC3DlL,IC3D2L,IC3D3L,IClD1E,IC1D2E,IC1D3ii,xC2D1E,IC2D2E, 

2 IC2D3E.IC3D1E, IC3D2R, IC3D3E,ICOL 

FORMAT  (18  rl'^rAE, II, A3, 13, 13,  18(11)  ,2X,I  ) 

GO  TO  10 
END 
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APPENDIX  5 

PROGRAMS  EDIT2  AND  EDITS 


PROGRAM  EDIT2 

LIS’;:  o :u>  , rr  (D:  4) 

^riEHSIDS  I0ffE(D:4)  , I HP  (0:4)  ,I«TN(0:4),TTS(0:4)  , 

2 in?I(0:'4),IOFF{0:'4),mDl  0:4)  ,ir1O2(0:4)  ,I:1d3(D;4)  , 

2 ir2DM0:4)  ,rC2D2(0:4)  ,10  20  3(0:4)  ,IC3D1(0:4)  ,103  02(0:4)  , 

2 T0303  (0  :4  ) , TOOL  (■»  :4)  , TOODSr  (D  :4) 

0> IL  IFIIE  (21 , ’STP'T. ' ) 

0!)LL  OFILF(22,  ') 

1 i = ) 

Of  21  , 'i,''SO  = HOO)  ID'  OF  (1)  ,IHP  (M)  , If!IN(fl)  , LINE  ( Jl)  , I TS  ( 3 ) , 

2 'T  (H)  , TON  (Ml  , lOFF  (3)  , I01D1  (HI  ,10102  (H)  , I0l0  3 (H)  ,10201  (3)  , 

2 r02->2  (3)  , 1020:’  (3)  ,70301  (H)  ,10302  (Ml  ,1030?  (3)  , lOOL  {(!) 
F0O3'»r(I?,'>(I?),'2,I1,?3,2(I3)  ,’4(I1),2X,n> 

2 00203=2,4 

RE*  0(21 ,5,ENO  = 23)  LOME  (3)  ,IHR(3)  ,I3IN  (3)  ,LTNE  (3)  ,TTS  (SI  , 

2 :T(3),IOV(*I,IOFF(3),I0  1d1(3),I01D2(3),I0103(3),I:2O1(3), 

2 10202(3), 10203(3)  ,103  01  (3)  ,10302  (3)  ,10  3 03(31  ,I0DL(3) 

!f=3 

IP  ( ( (I3R  (HI  *S0  + r3IN  (3)  ) ” (IHR  (3  - 1)  *604.imih  (3-1)  I ) .OT.  2) 

2 30  ro  25 

TP  (I ’■N'"  (3)  . NT.  T THE  (3-1)  ) 30  TO  25 
IP (~r  (3)  . NE.  (3’1 ) ) 30  "0  25 
IF  (IT?  (3)  .N’:.’’!'’  (3-1||  30  "0  25 
IF  (i:OL  (3)  . ’'0.  I 'OL  (3- 1)  ) 30  TO  25 
IF  (3.L".'>)  30  '0  ■'0 

’•P  (I 'OL  (3)  . ^0.  I^OL  (H-2)  . OP.I  OOL  (3i  . EQ.  lOOL  (3-3)  ) 30  TO  25 
"’0  '■ONTIN”E 

30  -^0  24 
->3  IF"  = 1 

24  IP  (K.  EO.  4)  ^ =-: 

25  00  3'!  3=1  ,<  -I 

I:0'7S’’(3)=:'101  (3)*100000000+I01D2  ( 3)*10000000tI''1D3  (3)*1  000000»- 
2 I'~?0''  (3)  T'n-'‘'T>4-ro202  (3)  *13000  + r02O?(3)  *"'000  + I03o1  (3)  •lOO^- 
2 10’02  (3)  T 1OH0302  (3) 

I?(I'OON*'(3l.'’0.':rO'TS-(3-‘').?.ND.ION(M).EO.TON(3-1))29,25 
25  Tf  (H.  l:  . 7)  -o  73 

:F(T'O'TNT(3).^0.-0OgNr(3-2).BND.I0N(H).E2.ION(3-2))29,27 
•>7  I’’(I'0nS'(3l  .'’0.I“O'IN"(3-3).»ND.ION(3).E3.IOS(3-3))29,23 

">4  ioohn:  (3)  =2 

TOI”*  (3)  =0 
I'1D2  (3)  =0 
i:*'  D?  (3)  =0 
10201  (3)  =2 
I'OO?  (3)  =0 
10  20’  (3)  =0 

T-7T1  (■*)  =2 
•^73n2  (3) 
lO’D’ (31  =^ 

'OH  (3)  ="* 
lOF” (3) =0 

T0  ITO’:^'  I.=  ’”‘00’.  I>I'0O7"  (3) 

0OOP  = ''O'^?+-:01^1  (31  '■101  ”2(3)  '■10103(3)  -Hoi'll  (3)  '•r0202  (3)*- 
2 1020?  (i)'-IO?Dl  ("I  ♦:-?02(3l  '•10303  (3) 

*'ON='Ot>ION  (’*1 
T0F”'  = rOF’’+’’0'’F  (3) 

?0  'ONT'S'TE 

rF('"-  (T)  .PO.'^"''’.’’SO.I''"  (1)  .EO,  1)  40,50 
IT  I p (3  ' O'*  L.  SO.  1 000000  0)  '90,150 
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PROGRAM  EDIT2  (CONT.) 


60,  ■’0 

6"t  !'■  (rro:"5 L . -’3. 1 n 1 1 looo*  ioo,i60 

■^0  IP(?r(H.F0.'P'"'.iN0.TT5(1).EC).‘’»B0,90 

po  i9o,i?o 

'’o  ipr'’r').F;5.»LPV'.!tn.TT'^(i),E2.i>  no, in 

m iPfiroriL.^'o.  imnn*'*  no,io'> 

in  IFf"r(1>.r:>.'LFVi.?v:».TTS(1).FQ.2)120,600 

.;>■)  TpfT-^-^L.Eo.nmnDD)  no, no 

no  I 0FF  = I M”  f rOFF 'HOOF  *2«-.  6) 

fON  = ’’Nr(rcK/''''or*2^.T) 

■JO  •■>  no 

15-!  I 0=’F=rN’T  fTOFF/f'P'^F  *'4  + . S) 

L0’4  = inrr^oK/nno?*'4  4-.5( 

10  "0  no 

ni  ioFF=ifi’’ rroF’'/poop*5«-.  F) 

L0M  = r”'^  fT’ON’/DOO'’*5«-.  51 
10  “0  '00 

n't  L0PF='^f:r  (?OFF/POOF*’f.  F) 

LON  = '^'!T(TOn/mOF*3  + .5l 
10  ’"0  ■’"'0 
no  r oFF=T’''p? 

LON=~OV 
00"  ’’K=F-‘' 

W!>I-  R {’2,Fl  ID’-  "E  (1  I ,1.1-  (n  ,rniN  (II  ,LIHE(1»  ,123  (1)  ,0’-  (1l  , 

■>  Lnl,LOFP,:r 

1°'  (IFI.  EO,  1)10  ro  FO'^ 

3-^  s L = 0 
■>0"'F='' 

1 0 fi  = '- 
”OFF='' 

IF  (K. so. 5l  10  -0  1 
’’O'  r^’fii  =i'"''’'i(F) 

IF"  (II  =I1F.  (K1 
rii'rc)  =i.irF(‘<’) 

LI*I-  (11  ^LIN’"(^1 
ITF  C 1 =ITS  (FI 
i:’('i  =-”(Ki 
lOM  (1  1 =ION  (K| 
lOFFC  1 =IOFF  CO 
I^ini  f =rini  (K) 

i:n?  (11  =r-ioo  (p) 

r?i D? (II =r-i D?  (Ki 
ri2i^i  Cl  =ii2n  (Ki 
r?2''2  (11  =r:2no  (P) 

11203  (1 1 =1:20^  (iO 
:73m  (n  =i-:pi  (K) 
r:3n2  (11 =11302  (Kl 
II3D3  ("  ) =i:?D3  (F) 

I'Di  (1 1 =rzoL  (?:) 

10  ”0  ? 

0 T E 1 0 
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PROGRAM  EDIT3 


Program  EDIT3  is  similar  to  EDIT2  but  is  used 
for  trains  opening  doors  on  both  sides. 

rJT"Gr=  LI\>E{.0:^),CT(0:43 

T T 'ir  VC  i:>.M  I -lA  TF  ( 0 :4)  , I HR  ( 0 : <»  ) , I M I ( 0 : ^ , I T S O : A I , 

? I'',N  (^:  A ) , lOFF  ( 0:a  ),  K.  1 niL  ( A)  , IC  1D2L(  n;  A)  , ir  103U  OtAi , 

’ TC  2^11.  ( 0:A)  ,IC202L(  OrA)  ,ir2n3L(0;A)  ,IC3'^IU0:A),  IC3D2L(0:A), 

2 1 ^3H3L  ( 'l:  A ),  ICiniR  { o:  A ) , 1C  10  2'’  ( 0:  A)  , IC  10  35  ( 0;  A)  , IC2D1R(  0 :A>  , 

5 ir?O?5(T:A),IC2H3P(n:A),IC3DlR(0:A), 1C3D2R ( 0 : A ) , I C3  03R ( 0 : A ) , 

■>  1 ROl.  (0:  A ),  ICC'HNTn  : A)  ' - 

riLLlFR  '^(21,'STRAI.M 
-M  I '-'^IL'"(22,  'STRA2.  •) 

'1=1  

( 21 , 5,  END  = F 00)  I?/'  T=(  •')  , 1 HR(  ,I  MI  M M)  , LI  NE  (M)  , ITS  (V  ) , 

5 '■T  ( v|,  I rv(  y)  , iPFF  (V,  , , icioiLIM),  IC10  2L(M  ),  IC103L(  M)  , IC20  1L(  M)  , 
2 K 2"^?L  ( M)  , ir  20  3L  CD  , I C 301L  (M)  t I C3  02  L ( M ) , I C 3 03  L ( M ) , I C 1 01  R ( M ) , 

? T*"  1 02h  ( M ),  IC  lO-’RI  -»  ) , IC  20  F ( ^- ) , IC  20  25  I M)  , IC  2035  (M  ) , IC  301P  ( M)  , 

2 IC  ?0?5  ( M ) , ir,3T3R  { M)  , If  CL  (V  ) 

POP'MT  ( IP  , ->{  12  ) ,a2,  I 1 ,A3,2(  I 3)  , 9(11)  .on  l)'72x\  fl) 
on  20  ••  = 2,  A 

RR-'  0 (2  1 , 5 , E0n=  23  ) lOATT  ( A j , IHO  ( N')  , I MI  N(  M)  .LINE  ( M)  ,l  TS(  M)  , 

5 CT(  M)  , I PM  , I ppr  (M)  , 1 cini  L(M) , IC102L  (M  ),  IC103LI  M),  IC2D1L(M), 

5 1C202L  (M  ),  IC203L  (M  ) , IC  30  IL  (M  ) , IC302L(M)  t_IC  303L(M)  , _ 

2 ir,  lOiri  ( -j<)  , I Cl  025  CR  , IC103P  (M) , IC201P  (M  ) , 1C2D2P  (M  ) , 1C2D3R  { M ), 
p jn  ?nip  ( -R  , IC302R  (M)  , If  3033  ( M)  , IC0L(  M) 

K = 

1 5 ( ( ( IHF  ( '1)  *A0+I  («  ) ) - ( I HR  ( M-1  ) ^60+1  *^1  N(  M-l  ) ) ) .GT.2  ) 

? C.n  TC  2 5 

ic  ( L INF(  M ) . L I NE( 1 I ) G("  TO  25” " ^ 

IF(rT(M)  .NE.CTCM- 1)  )C-n  TO  25 
I5(ITS('").NE.ITSC^-1))G0  TO  25 
1=  (I  COL  (‘1  ) . EQ  . ICOL  ( M- 1)  )G0  TO  2 5 
15  ( M.i  t.3  )G0  to  20 

I5(  ICPL  ('^  ) .EC.  ICCL  ( M-  2)  .OR.  1 COL  (ii)  .E_Q.  IiIQ.LiM-3J  1 GO  TO  2 5 

Cnr’TIMlE 

GO  in  •JA 
I R I = 1 

I F( K .FC. A) K=5 

OR'  3'^.  ....  ... 

1 rnUNTf  VI  ):=  1C  10  IL  ( “)*  lO**  17-H  (:Id2L  ( M)  *lb**16+  IC1D3L  (M)  *10<=*  15  + 

2 10  2D1L  ( M )*10  + * lA  + IC  202L  (M  ) «10**  13+IC  20  3L(  «)♦  10**12+ 

2 1030  1L(M)«io<:«ii-Mr3  02L(  M)  * 10**  10+ 1 C303L  (>'  )*10**9  + IC  IDIR(M)* 

2 1 0**8+ IC102R  (M  ) *10**7+IC  10  3R  (M  ) *10**6+ IC  20  IR  (M)  *10**5+_ 

2 I02O2R(  M)*lo**A+ir.203P  (M)*10**3+IC3D1R(M  )*io**2+ 

3 JC302R  (_Mi*10+IC30  3R  (MJ 

I F ( I CCUNT( M) . EO. ICDUNT (M-1 ) . AND . I 0N(  M ) . E Q. I ON ( M- 1 ) ) 29 , 26 

I5(M.LT.3)G0  to  28  

I F ( I COUNT ( M) .EQ. ICGUNT (M-2 ) . AND . ION( M ) . EO . ION ( M- 2)  I 29, 27 

IP  ( I CO'UNT(  M)  .EO.  ICnUNT  (M-3  ) . AND.  I C N ( M ) . EQ.  ION  (M-_3  H 29 , 28__ 

1 COUNT! Ml =0 

J C101L(M)=0 

IC102L (M )=0 

ICl,D3U.M)j^ 

IC2blL(M)=0 

IC202UM)  =0  

IC2O3L(M)=0 

JX3.niJ.lMJL50 ^ 

IC302L (M) =0 

jr3  01LiMJ  = 0 . 

icioiP(M)=o 

IRl  Q2P.  IM  LpJI 
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PROGRAM  EDITS  fCONT.) 


- iciA^^-u)=a 

I r^niP  ('I  ) = n 

. LC2QZEI*iJ^-. 

IC?03P(«»=^ 

. IC^ni'^(Mj=.o...  ..  . : 

I'*  312P  (^  ) =0 

. . iC3n3°<M)  = o 

In^l(M>  = o 

IQF.EIMI-D 

1 TnTAL=I  T0T4L+ICrUNT  ( »') 

TnnR  = ncc)R  + iciDlL  {'*) +ICIDPUM  ).t.IQJJ)3L  LU >IC2DlLIMJ  + LC2I!2Ll!llt___ 
? IC.2D3l<'^l-HC3niL('^)  + IC302L  (*<)■►  I C3D3L  ('i ) + IC I DIP  (M  ) ■»■ 

7 K n2R(  M)  + IC133R('11  + IC201.P  (M)*I{L2.02Ki'll+K2D3B(^LLlL 

2 ir^Pl  r ( M)+ir.?D?P  (M) +IC3D3P  (M  ) 


''n^)=T^^'+I  0N(M) 
TTPc  =TCFF  + IOF_FlMi 
30  r'l^'TjMijc 


IF(^Tn)oFP.'PCC*  .ANO.ITS(  U ,E0.1)40, 50 
. . i=j_LT^TAi..=a.u i<'ia**ie>+u»iO’*L<=  1111.934,15.0 


5n  l = ( r Tn  ) . EQ.' PCC*  .AND.  ITSi  1 ) .EQ.2160,1.0  

60  I >"{  TTf.Tn  .EO.  ( I I’i-  10**  16+  11*10**  13+11*10**7  + 1 1*10**  A > 1190,150 


70  IPCTm.EQ.'PCC'.aN0.ITS(ll.EC.3)80,90 

80.  . .JP.ilTCTAL.  EQ.l.ll*10**l6,+  li,*10**13<^... 

r.  ll#ir*<rio+l  1*10**7+11*10**4  + 1101)  19  0,  170 


go 
1 m 

no 

120 

150 


1 60 


170 


1 0 0 


1 9 0 
30  0 


IF{rT(l).FQ.'LPV'.AND.ITS(ll.EQ. 1)100, no 

nnTCTAL.PQ. in  000000)  190,130  ...._.  

lu  + io**i‘5’+i'i'*'0'*+lp 

I F.(0T(  n . Eq..'LRVnAN3  .ITS(  11  .EQ42)  120,  530  _ 

n(  ITPTAL  .PO.  ( 11 1 111*  10**11+111111*10**3)  >19  0,170 
L^FPr  PIT  < TnPF/0GOR*4+ .5  ) 

L PM  = I \T(  T PN/D'^^R*4+.  5 ) 

00  rp  pop 

L^f^‘"  = INT(  TnPF/0Of^'^*0+ .5  ) 

l_p,i=  I Mjj  T0M/nnpp>^8+,  5 ) 

p n TO  T 0 0 

L'"F'-  = in^T(  TnFF/oooP*12+.9) 

I M(  Tr'j/opo-j*i2+.  5) 

p n T p TOO 

i.nrF=  IM  ( TOFF/nOUR  *6+ .5  ) 

LGM=n.T' TGN/D2CR*6  + .5  ) . - 

Tp  ppp 

1 nfi-TP^- 

I K = K-1 

vt>lTr(22,5)IDATP(l),IHR(l),I'^lM(l),LINF(l),ITS(l),CT(l), 

.3  LCN.Ll'FF , IK  . . ....,  

IF(  !F  l.FQ.l  )GO  TP  5^0 
npTAL  = 0 

nnr 

TP':=o 

rp  F F = o 

IFr<.FC.5)GP  TP  1 

IOA-^r(i  ip^TC(K  ) 

IH'’  ( 1 ) =IHR1.<) 
r*  I*'  ( 1)  = I Mr<(  K ) 

L ( 1 )^LI  \E  (K) 
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PROGRAM  EDIT3  (CONT.) 


iTC.  (1 

''Til  ) =CT  ( K ) 

!■"'■(  1 ) = 

T-irr(i  )=in:rp(^) 

T r 1 1 1.  ( n = I r 1 '>  1 L ( K) 

1 r 1 1’L  M ) = IC  lO’L  (*<  ) 

I- in-^L( n = icn? I ( K 1 

I'‘^'^IL(1  ) = I''’n]L(K) 
ir7-)2L  (1)  = ir21'>  L(  K) 
r'“'^3L(  n = K?''^L(K) 
I r ■>  "*  1 L ( 1 ) = I C 3 T 1 L ( K ) 

ir  ?L  ( n = ir3~)?L  ( k ) 
I r.->-'3',  ( 1 ) = ic-<T’L(  < ) 
ir  I"  ( 1 ) = in  3ip{  K) 

m (I  ) = icn23(K  ) 

I 2n3'^  ( l ) = ICl  ( ■< ) 

n p"'  IP  ( i ) = I c?mR  ( ) 
I r->^2P  ( 1 ) = I C?T2"  (K  ) 
ir?33'’(i  ) = ir?n^p  (ki 
I C 3"'  iK  1 1 ) = IP3~ilP  ( K ) 
n3^.?-’  (1 ) = I rn?f’.  ( K ) 
I r 3"'  3=  ( 1 ) = I''3"'3P  ( K ) 
I rni  ( 1 )=  ICP'L  (K  ) 

? c 

- \l-> 
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APPENDIX  6 
PROGRAM  SEPAR 


Program  SEPAR  sorts  the  data  into  two  inbound  and  outbound  files. 


D ? N F M • L . 

=■  '-IT. 

1 PaKLEC- 

■ipi 

Li:  T..1. 

-r'  p")nr  [_T  L?. 

31")  L‘:  . 

12%  ~H’K3E  HF'i'F  PTJP+12. 

31  -'i  L’. 

T •>  . p r P ft  •)  , 

F I D LIKE  ? c>  " L " 1 ~ ” 3 ~ " 1 r '"J  ” IE  " P " . 

3Y  :i01F  r'fM. 

’3TX)-  -)N  '1'TD’T  nil’  M7K  LIRF  Tr’INSIZE  ’Y?El'’''''r  C 

P"?Si:!  ?’-5'lP'='  ''OI'Frr  F'IT  lU  "^2  I2  n "3  I3  13  lx  II  S'!!. 
PT.’n  I • 31. 

31~"'  "Y  TO’’” 

TM'T  1 Pakriv. 

Ofry^-  ->  y'YD'Y  ill'  :il’1  LI)JF  "'AIN3IZF  ^Y?F1F'?.H  5 

ClIl^lT  FI  FI  AT  D'4  12  12  ^2  S 
■"1  ’’  13  13  iy  II  FSl. 

1 r"  *■ . 


6-H6-2 


.W'. 

-.v.*''r  vr  '•r:v;4;^.!t:i  ^r' 


V'  . ,5...  I 


Sstir'  fc 


.... 

• <W'%> 

'/•"  ■■  ' ’'jaar^ 


■ •'  .V4'  « ■'  f,'.  •• 
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FINAL  OUTPUTS:  PAKRIV  AND  PAKLEC 

PAKRIV 


(U 

4J 

CO 

0)  0) 

B a 

mH  *H 

0) 

cq  N 
Uj  -h 

c 

44 

44 

0 

a 

H J 

H 05 

0 

0 

1Q7Q052'* 

7 0° 

2L?7 

*5 

2 

2' 

19790^22 

710P 

21  RV 

14 

1 

2 

1079-1522 

7151 

•'L°7 

3 

0 

2 

19‘'9''S22 

71737 

2P77 

(3 

3 

2 

1 3 •»  a 7 52  7 

71877 

? P '*? 

3 

9 

2 

19"’B9522 

720F 

21  FV 

5“ 

1 1 

2 

^3792522 

722H 

1 P'7 

0 

9 

2 

197P3522 

7231 

1 P77 

3 

9 

2 

^■’30^22 

72537 

7^77 

0 

9 

2 

197B''‘i''2 

725P 

21  FV 

3 

7 

2 

19799522 

7 2 51 

1 P77 

3 

9 

2 

19’92522 

7261 

11  FV 

D 

3 

1 

13TP7577 

•'261 

1 3-- 

9 

3 

1 

197B"*522 

7261 

1 I 2V 

9 

9 

1 

19799522 

171? 

2L=»V 

C 

7 

2 

19’’'’''‘i22 

728'7 

2?'7 

9 

3 

2 

“'?''9a5?2 

7?:’1 

1 P77 

9 

9 

2 

I9797522 

73417 

’P'7 

3 

-) 

1979'»S'’2 

7 2 5 - 7 

9 

3 

0 

13797522 

7 4 0 2 

21 FV1 42 

34 

4 

13-.P3522 

74  1'7 

TL'^V 

9 

9 

4 

•19793522 

7481 

1 P'7 

3 

7 

4 

•!  3 -I  3 3 52  2 

740-- 

2^77 

9 

9 

4 

13797522 

7501 

1 P77 

4 

'97P0522 

75102 

•'  L’V 

9 

a 

■> 

1 97«3522 

7511 

1 T "V 

-5 

-) 

•'521 

1 P - 7 

9 

3 

1 

7521 

1 1 RV 

P 

2 

1 0 0 9 

•’52  = 

2 L ? V 1 1 9 

24 

4 

19733522 

7 531 

1 T.  ?V 

1 

1 0f  0 3 537 

7551 

9 

-3 

^9733522 

•'55'7 

1 J 
< 

7 

•% 

*> 

''56'7 

1 L’V 

9 

9 

1 

39-03572 

7 5677 

3T  cy 

9 

-> 

1 

19703520 

7561 

1 0 7 

3 

4 

33703572 

756  F 

21  '"V 

-) 

1 

3 

■»  97  a 3 52  7 

7 C 7E 

•>L”V 

-9 

9 

7 

33700577 

7511 

1 ?" 

2 

i07C'^5*'-) 

1 0'7 

>2  n 7 - 

9 

3 

4 

19‘^3''B22 

1 1 

’I  ?v 

■5  *3 

2 1 

4 

1970757'' 

9 41 

1 P" 

3 

9 

4 

13--03522 

1 5'7 

' 0 P" 

9 

*> 

u 

1 D T 2 5 0 

P 517 

■3077 

j 

9 

4 

39793577 

1 61 

7107 

•> 

3 

1 

-10-p'^n;-)') 

1 61 

■'  L3V 

9 

9 

7 

3-»or>q^p 

9 71 

1 I ?V 

*> 

1 

1 9-’035-f 

a c ■< 

1 L’V 

9 

0 

4 

13-'9n527 

9 •(  n 3 

■'  P'7 

3 

2 

l?‘'  = 7‘^22 

p T 0 a 2 

l_2_7_7_7 

3 

9 

'1 

^ :)-TQr\c-)-) 

112'' 

2 T f V 1 2 1 

4 

° 1 3 77 

' 7 - 

9 

4 

1979''5'I7 

1141 

1 ' ?V 

3 

1 

1?'’50527 

9151 

D 7 7 

3 

3 

■> 

39-.q-52f 

31  5~7'’p-~ 

9 

'I 

1 

•'970752'' 

11  51 

•1  -3  •^  — 

9 

3 

2 

1370-1522 

1161 

1 P77 

1 

^ ra-f  0 ^ c - ^ 

8' 61 

in'*- 

3 

3 

*2 

19‘’B''522 

9 3 fr 

•'I  ''V 

3 1 

4 1 

u 

•'2-0352'' 

9073 

« n - ^ 

3 

9 

4 

124P 

3t  0V3  30 

35 

4 

1 c ^ 0 ^ S ■>  ? 

0 -ir  1 

1 ?77 

•^9 

r\ 

4 
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Surveyor 


PAKLEC 


I 


Date 

Time 

Line 

Train 

size 

c 

0 

Off 

If  of 

1 q ->  5 7 'i 

730-;- 

IT.  7V 

9 

2 

13'’9352'4 

72  1 M 

21  = 7 

-5 

35 

“> 

1979''S2U 

73592 

1 P-? 

0 

0 

2 

73P0 

1 L77 

0 

0 

1 

19Tq«,52’4 

73OL 

20-21 371£2 

2 

1 97  0-I5T-J 

7ir7qr 

0 

0 

1 

1 97  9 ’>S2'4 

741 1: 

1 1 R7 

-1 

0 

2 

1 9 ^ B 2 S ■>  ’4 

742'1 

2LOV 

75 

73 

2 

1 q 7 p ''  5 2 '4 

7 4 31 

3P-? 

2 2 

43 

2 

197B2B24 

7449? 

^ ® rr 

0 

0 

2 

197p7<i2-4 

74  PM 

2L  B.7 

24133 

1 

10799^714 

748L 

2L7V 

80758 

2 

1979957-4 

74P9? 

1 P?? 

0 

0 

2 

ia-905'>u 

7499? 

-ID-.- 

0 

0 

1 

19’'R952'4 

7 4 99? 

' P?? 

q 

0 

1 

75OL 

■ L7V 

5 f 

5 4 

1 

19797524 

753N 

2T  P71 

''5213 

2 

99790524 

755L 

2 0 2 

7 5 

57 

2 

1 q-'9057!4 

7 5 7M 

21 F7 

151''1 

2 

19799574 

75P9? 

• ^ M 

0 

0 

2 

1979757-4 

75BM 

2 1 ° 7 

1 

26 

2 

197095711 

3 ON 

r !?  y 

c * 

7C 

3 

1979757-4 

3 4 9? 

' P*  - 

1 

0 

4 

3 ->  0 7 5 7 14 

3 4L 

-?'? 

437  0° 

3 

1 9 ' 3 0 5 2 4 

3 51 

op-- 

15 

6 5 

4 

7Q-'  = q57f4 

3 ■^■9-' 

L?7 

0 

0 

4 

19797574 

3 fif 

?P?- 

3 

54 

4 

1979757-, 

3 79? 

. D - ^ 

a 

0 

4 

1 9 7 q 5 2 4 

3 R9- 

'*  n-*  « 

4 

^97 9 05“ 4 

3744 

- r_  0 V 

571 

P C 

4 

19797574 

3151 

7p-- 

5 

QO 

4 

1979757-4 

3-6L 

2P?? 

4 5 

78 

4 

1979757(4 

P179? 

1 7?? 

3 

1 9 7 Q Q c 4 

81BL 

D D 

0110 

4 

1 9 T 9 7 5 ■>  -4 

31  9N 

1 r 77 

2 

0 ^ 

4 

'9"'905'4 

3 2 0 9? 

- 0^2 

0 

0 

4 

197  9 ->5^  4 

3249- 

l_TLZ_r 

»> 

4 

49730524 

2p-? 

27 

4 

1 9 7 q ■>  5 2 4 

327' 

7p?- 

1 

5 2 

4 

19707C7I^ 

3 2 9 - 

• 71  — 

0 

*> 

4 

19737574 

329N 

1 r ?7 

1 1 1 1 

4 

19797574 

8 7 0 9? 

1 D ^ 

0 

a 

4 

19737574 

3 ’ 7 D 

1 r sy 

0 

1 

19’075'4 

3775- 

- p - 

0 

4 

19-37524 

33  29? 

2 p ■'  r 

3 

4 

19737c- 4 

3'5L 

D r>  — 

0 1 

4 

19597524 

3 3 6 N 

7 T P 7 

217  7 

4 

19 ■’3  0 5?  4 

3 G 9 ? 

0 

0 

4 

19-97574 

378' 

2P-? 

34112 

-a 

7 4 0 3- 

^ n -» -* 

0 

9 

4 

3 4 ' 

■5  D - " 

4 

95 

fi 

<3-o-7t;->i. 

R U ? 'I 

' r.77 

0 

3r 

4 

1 9 - 0 7 5 7 ^ 

B 4 .7  9 - 

ip-- 

*) 

C tr.  4 

9 '4  f r 

r>  — — 

q-4  1 

1 0 

4 

1 Q -»  0 ^ C ^ J 1 

94-p 

1 r - V 

-1 

1 

94'’9? 

- • D — 

0 

■3 

1 3 - 3 7 C - -J 

34  3 N 

1 » D y 

*3 

ar. 

4 

1 ^ D 0 ">  U 

757[ 

0 

5' 

a 

^q7r>'^c-*4 

3529- 

ip-? 

*s 

1 

ti 

15-3752-1 

?77L 

- D — 

1 - ’ 

24 

4 

1 4 -r  n ->  4 

3529? 

2 f’  ? ? 

■3 

CO 

u 

o 

>% 

Q) 

> 

u 

3 

cn 
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APPENDIX  8 
PROGRAM  PROPFL 


Program  PROPFL  generates  the  actual 
data  for  individual  stations. 


time  series  of  the  observed 


d:'!Fn-Ion’  :h'-  tsO'’')  , i :;:m  (sor)  ,icF(5f;0)  ,iorF(50C) 

:-v)'jtLE  r*EviiiZo:;  fileo'j,  fil li:;,  di:,ct?v: 

7  7 F 7 5 

5 1 u: : r ;f7U7  file  e -•,?) 

ACCli^T  o,  FI  I FIN 

TYPF  4 

4 FOI  EM  ( ' GL'IPFT  Flli]  N’ilE  = ' , S) 

ACCEPT  ■j.FiLFOU 

6 rCr'lAT{AlO) 

urn.  (U:iIT=;i  ,DFVICF=  'uFF'  , ^C'~F,FS='  FSGIN’  , FILE=FILEI:i) 

OPiN  rJNI  T = 21  , DEVIC^=  ' FrU' ' , 7-CCISS=  ' F ZuCUT*  , FIL '■=  FI LZO’J) 

TYPF  7 

7 F0F.1?.Z(’^0  Wh.M  D_>  ZCTION7  ' , IX  , ?) 

ACCEPT  6, JIT 

8 FC  RE  AT  (AID) 

TYPF  r 

9 FCr..AT  ('AT  WHAT  S i A T Z O N ? ' , 1 X , S) 

ICCZPT  IO.FTFT 

10  rCF  f.  'T  (A10) 

wr.lTZ  (22,  11)  STA'^MIT 
h-ITE  (22,12) 

12  FOR  FA-  (//2X  , IF.Z'  ,2X  , ’ HEA  D.;  A Y ’ ,2X,  ’ FAFS  CN",1X,'QN  FATE’, 

2 IX, 'PASS  OFF', lx, 'OFF  1..7TL  ' , 2X  , ' ' RAIN  IN  bZT') 

N’  = '' 

i: ZEr=0 

DO  ICO  I=1,5t'0 

ZLAD(21 ,20. 2ND=T00)  IllT  (T)  ,IniN  (-)  ,T5  AIF  (I)  ,10F  (I)  , 

2 lOrrME) 

20  FCr.EAI  (8X,12,I2,  '2  ,'4X,I3,I  j,  2a) 

?AS5o:;=:oN  (_) 

PAS0FF=I0FF  (1) 

NTJi-;=Z-N-1 

IF  ( DIF.  CO.  ' LECHEFF  £' ) SO  10  23 
IF  (TI  A.IN  (I)  . ZO.  ' F.  ')23,1J0 
23  IlEOWA.  Y=  (IHr  (Z)*60  + IE:'i(:))-ZTIF.F 

PEDWAY=HEDJA  Y 
IF  (hlDWAY.FO.O)  :>EDWAi=1 
?A-FGN  = P.f  SSO:; /PEDWM 
I'.ATOF  F=?ASOFr/?_  FWAY 

Wr.ITE  (22 ,50)  ZHF  (Z)  , liHF  (I)  ,iiEDW  AY  , ION  (I)  , F.MEGN.IUFF  (I)  , 

2 PAZUFF.NU:'! 


11 

FOr  ilAT  ( 1 X,  ' AT  STATTON=' 

, IX,  A 1C, 'GOING  TO',  IX, A 10) 

50 

FORMAT  (2X,12 ,Z2 , 4X , F5.  1 

25X.Z3  ,5X,F4.  1.8A.Z2) 

60 

ITIilS=ZtIR  (I)  ^60  + IMZ:i(I) 

S = I 

100 

CONTZN’DE 

500 

END 

8-1/8-2 


#^--'siiP,>.' 


" " ■v'/;:iJi 


i;!f: 


■ ■ i : i kf 

''  4'ii  i 


t. 


■a* 


m'i' 


# 


. _.  .._.  , 


pr' 


. . .1  -V  ^ 


ii''v.  ,*'v^  ..  .*  r ■» 

.1  •’•  • i.  ..1  .», 't..'«,:is  <v  ,11.'' 


''■■  ',:'r  '■^, Vi.  «5i 


■■  ' 4 


■ .f f'v  4|. M-j-iP,  :■••. : fc  «-5^S»f|»^(^  riiM 


■,S  ..  ^ : iilki  ‘ 

% . f '’f  '.  "X;  ,,.'  •"»; 


>>2'' 


m 


<U  ^v, 


»-  * '■'  ..'HI  »*:¥,? 


i'll: 

Ivifr  ■'■C' 


',#'  M" 


■«  '■»»»  11 


APPENDIX  9 

STATION  PROFILE:  PAKPFL 


AT 

STATION 

= PARK 

ST. 

GOING 

TO  RIVERSIDE 

'Mr 

5 E D 4 > y 

P'a'  014 

0'1 

2 S.  r E pas 

S OFF 

OFF  RME 

r 7 ».  I N I H 

7 n 

4:0.0* 

3 

0.  !0 

2 

9.0 

0 

710 

10.0 

' 4 

1 . 4 

a 

0.  8 

0 

72? 

10.0 

3 3 

3 . B 

1 1 

1.1 

3 

3.  0 

9 

1 . 8 

7 

1.4 

a 

7?7 

2.  0 

5 

3.2 

7 

3.5 

4 

7!i0 

M.  0 

1 4 2 

10.9 

34 

2.  5 

4 

7 '5 

17.2 

1 1 0 

8.5 

24 

1.8 

3 

0 

2 

0.7 

1 

0.  3 

5 

7M 

1 . " 

2 

2.0 

0 

0.  0 

0 

0 1 

4.  0 

a a 

7.  ? 

31 

7.  3 

2 

a 1 2 

11.0 

122 

10.9 

21 

1.  9 

9 

a 1 

3.  0 

0 1 

19.  2 

41 

8.  2 

7 

a-7  U 

1 . 0 

118 

15.  9 

35 

5.  0 

1 

5 '4  0 

1 0 

04 

3.  9 

33 

1.9 

7 

a '4  7 

7.  0 

40 

5.7 

42 

6.0 

9 

a *j  a 

1 . 0 

3? 

77.0 

21 

21.0 

0 

054 

5.0 

1 5 

2.7 

3 

0.3 

4 

3 0 

15.  0 

?.  6 

1 1 

0.  7 

20 

?1  1 

7.  0 

25 

12.  5 

2 

1,0 

3 

- , c; 

4 . 0 

a c 

7.  8 

8 

2.  0 

1 

a 2 Q 

17.0 

7 4 

1 . 8 

3 9 

3.  0 

9 

: ~4 

5.  0 

■^  8 

5.  3 

23 

3.  9 

5 

a '4  2 

5.  0 

a 1 

7.  5 

6 

1.0 

5 

10.  0 

7; 

7.2 

57 

5.  7 

7 

5.  0 

17 

2.  8 

5 

0.8 

7 

10  12 

15.0 

- Q 

2.  4 

1 0 

0.  6 

f 

1 1 ">4 

■>2.  0 

7 4 

5.  2 

p 

0,0 

a 

1 0-^5 

1 . 0 

' 2 

12.  0 

c 

9.  0 

2 

1 ''M 

1 . 0 

1 2 

12.0 

14 

14.0 

7 

aatj 

188.0* 

^ -N 

0.  4 

72 

0.  1 

5 

■”4  2 

8.  ^ 

2.  5 

15 

1.9 

4 

■ ■'OS 

”.  0 

' 2 

5.  9 

79 

3.0 

2 

1 a r Q 

4 . 0 

01 

2.  3 

3-7 

3.  3 

2 

'45 

0 

1 5 

2.  5 

13 

3.  0 

4 

1412 

7.  0 

0? 

10.3 

24 

3.  4 

1 

14-2 

r . 0 

04 

2.  2 

6 

0.5 

3 

1 4?P 

5. 0 

21 

4.  2 

1 8 

3.  5 

5 

■"Oil 

1’.  0 

1 7a 

7.8 

21 

1.6 

7 

‘ 4 8 

7.  ^ 

05 

2.6 

4 

0.6 

4 

’4  5? 

0’  . 0 

a - 

* . 9 

a a 

3.  0 

6 

1 3 7 

4.  0 

a n 

a ^ a 

7 

1.8 

3 

"O'? 

- . 0 

1 ■'  7 

1 ? . 0 

78 

8.  7 

1 2 

13  7? 

1 ’ ! 0 

7 2 

3.  4 

68 

6.  2 

1 1 

1 5'’6 

^ a 

a r 

?’  . 7 

86 

12.  0 

2 

M ?2 

5.0 

1 1 

1 . B 

3 0 

5.  0 

4 

• -a4- 

T . 0 

' 4 

. 0 

4 

7.  0 

3 

10^7 

a 0 

16 

5.  ? 

14 

4.  7 

2 

^4C 

12.  0 

<;  a 

3.  7 

1 7 

1.4 

5 

■'305 

5.0 

4 5 

7.5 

33 

5.  5 

6 

15  5 

M.  0 

3.2 

2.  0 

15 

1.5 

6 

15  10 

10.  a 

3 0 

2 . 9 

2° 

2.  7 

8 

■>  r 1 a 

0 

4 P 

'5.0 

57 

1 0.  0 

0 

1 5?'' 

a . a 

5 0 

5.  7 

29 

3.  2 

p 

* 3 M 

a a 

a a 

10.'' 

21 

7.  0 

a 

154  1 

11.0 

22  1 

20.1 

1 8(2 

16.9 

a 

'632 

1 ' . 0 

100 

7.  ? 

1 Q 

1.7 

1 0 

1 5^5 

■ . 0 

5 5 

22.0 

14 

4.7 

3 

* large  headways  indicate  break  in  data 
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STATION  PROFILE:  PAKPFL  (CONT.) 


1 '11 

IP.  1 

111 

5.9 

0 

3.  3 

“f  71  P 

7.  3 

1 1 ? 

15.1 

51 

7.  3 

1 121 

3.  T 

1 3P 

52.7 

92 

30.  7 

1 ■'"'fi 

5.  3 

5 6 

13.  2 

1 7 

3.  '4 

1 13D 

'1.3 

2P 

5. 9 

16 

4.  3 

' 7’0 

3.  3 

29 

23.0 

12 

12.  3 

17^1 

1 . 7 

'4  1 

’4  1,3 

23 

29.  0 

173? 

7,  3 

72 

11.0 

40 

70.  3 

1 7’4  2 

3.7 

P7 

9.  7 

Q 

1.3 

1 7 U U 

7.  3 

0 

4.  5 

4 

2.  3 

1 p c 

15.7 

1 1 P 

7.  4 

3« 

2.4 

• P 6 

5.  3 

3 7 

5,  ? 

49 

8.  3 

1 p p 

?«  7 

? 3 

15.3 

2 

1.3 

^ 31  0 

7.  7 

ttP 

23.0 

7 

3.  5 
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STATION  PROFILE:  PAKPFL  (CONT.) 


AT  STATION  = PK  ST.  GOING  TO  NORTH  STATION 


E 

H~;  uv.i.  Y 

p ’ s 2 4 u 

.1  ~ - 2 E 

PASS  OFF 

OFE  F.ME  : 

IK  lU 

73  1 

451  . 

7 5 

: . 2 

36 

C.  1 

1 

743 

11 . 

•) 

2.3 

7 8 

7.  1 

5 

748 

6.  0 

2 4 

4.  0 

136 

2 3.0 

2 

753 

r,  r\ 

Mo 

21.2 

21b 

43.0 

5 

757 

4.  0 

1 5 

3.8 

10  1 

25.  3 

1 

75H 

1 . c 

1 

1.0 

26 

26. C 

1 

J 

2.  0 

53 

2 6.  5 

7 9 

3 9.5 

0 

n 1 4 

14. 

6? 

4.  5 

109 

”7’.  8 

7 

b1  i 

' . 0 

C . 4 

9 0 

1 8.  0 

4 

o29 

1 3 o 0 

1 

1 

1 1 1 

11.1 

5 

c2f, 

7 . 0 

2 

0 

1 3 0 

1 8.  b 

5 

f<4  3 

^ n 

ij. 

. c 

6 8 ' 

1 4,7 

4 

ir48 

(.  . 0 

3 

0.  5 

4d 

7.7 

4 

65o 

b.  ( 

1 2 

1 . 5 

b 9 

11.1 

-J 

5 19 

2 3.0 

1. 

3.  0 

26 

1.1 

14 

922 

9 . 0 

7 

2.  3 

4 6 

15.3 

5 

97  8 

i =0 

■>  3 

2.  2 

75 

12.  5 

6 

93  C 

2 . 0 

6 

31 0 

4 2 

21.0 

2 

9 40 

M.  0 

0 

0.  0 

54 

5.  9 

8 

94  9 

9 . ' ' 

6 

C . 7 

42 

4.7 

7 

95.4 

1 o 

70 

14.  0 

2 

95  9 

5.( 

1 3 

2.  6 

0 (■ 

16.0 

2 

1 0 9 

10.0 

0 

0.  0 

2 4 

3 « 8 

5 

1C  14 

5.2 

1 1 

to  e ^ 

46 

9.  2 

3 

1027 

13.  2 

M 

■I . 3 

59 

4.  5 

6 

U 39 

12.2 

6 

".7 

5 9 

4.  4 

5 

10  5 1 

12.  0 

c 

0.  8 

1 1 5 

9.6 

3 

1 C 5 B 

T A 

r 

1.  1 

4 9 

7. 

2 

1 524 

2 6 b . 0 

3 2 4 

1 . 2 

3 1 2 

1.2 

4 

1539 

f . 

3 2 

5.  0 

6 2 

1 . 2 

2 

1 5a3 

12.0 

: 4 

2.6 

4f 

3.5 

1 1 

1 544 

1 . ■' 

o 

6.  0 

7 9 

79.  0 

1 55  0 

6.0 

1 4 

2.3 

9 3 

15-  5 

4 

U 9 

•: 

C . 0 

27 

2.7 

M 

1610 

1 0.  J 

0 A 

5.  2 

20c 

20,  6 

1 

1o1  6 

; ^ i*' 

39 

4.  9 

3 

4.  9 

4 

16.1 

3.0 

A A 

7 ' 

74 

24.  7 

5 

1627 

>:  . 

14 

2.  3 

1 1 4 

19.8 

7 

16  3 4 

7 . J 

1 4 

2.  7 

1 3. 

1 8.  9 

6 

1o37 

J • ? 

1 1 

3.  7 

5 3 

17.  7 

2 

1646 

11.0 

? 

0.  7 

52 

4.7 

4 

17  0 

1 2 . 0 

2.  9 

1 0 3 

3.6 

7 

1713 

13.  J 

34 

6.  5 

1 32 

10.2 

1 O' 

1719 

4 5 

6. 7 

1 0 4 

17.3 

3 

1 72  5 

■..0 

1 7 

2.  8 

6 1 

10.  2 

4 

1727 

. . ■' 

1 3 

6.5 

3(. 

1S.0 

1 

17  4 3 

U.  0 

J 

0-3 

7 4 

4.6 

1 1 

1 744 

1 

1 2 

12.0 

56 

5b.  0 

0 

1o  9 

7 5.  0 

4 

0.  2 

62 

2-5 

1 1 

1 (5  2 8 

10.: 

45 

2. 4 

99 

5.  2 

9 

*Large 

headways 

indicate 

breaks 

in  data. 
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APPENDIX  10 
PROGRAM  CONTRAC 


Program 


5 

6 
7 


1C 

11 

13 

14 
12 

15 

16 
18 


2C 

21 


24 


25 

3C 


CONTRAC  contracts  PAKPFL  into  half-hour  intervals. 


DIMENSION  IT  IME  I ICC)  , INK  100)  , I CM  100)  , I OFF  (ICC), 

2 TCM50),TCFF(50),TOFWY(50),FECWAY(1CC), 

2 RACM5C),RACFF(50). 

DCUELE  PRECISION  FILECU,  FILEIN 
TYPE  5 

FORMAT! ' INPUT  FILE  NAME  = ',$) 

ACCEPT  6,  FILEIN 
FORMAT (AlC) 

TYPE  7 

FORMA T ( ‘OLTPCT  FILE  NAME  = •,$) 

ACCEPT  6,  FILEOU 

0PEN(LNIT=21,DEVICE=*DSK',ACCESS=*SEQIN*,F1LE=F1L£IN) 
OPEN (LN  IT  = 22,CEV  1CE='DSk  *,ACCESS=  'SEUCUT*  , F I L E = F I LE CL  ) 
DATA { INTI K)  ,K  = 1 ,28) /700  ,7  30  ,800 , 8 3 C , 9 00 , 93 0,  1000, 10  30  , 
2 1100,113C,12CC,1230,13CC,133C,14CC,143C,15CC,153C, 

2 16CC  ,16  3C  ,17CC  ,1730  ,1800,1830,1900 , 1930 , 2000 , 2C3C/ 
K=1 

FCRMAT(////,2X,14,4X,F5.1,5X,I3,14X,I5,19X) 

DC  25  T = l , ICO 
IF(  I .EC.l  )13,14 

REAC(21,1C)ITIME(I),HED^AY(I),ICMIJ,ICFF(1) 

GO  TO  15 

R EAC  (21 ,12  ,ENC  = 2C  )ITI  ME(T),HECV^AY(I),  ]CN(IJ,ICFF(IJ 
FCRMAT(2X,I4,AX,F5.1,5X,I3,14X,I3,19X) 

I F (I  TIME!  I ) .LE.  INT(K  ) )16,21 

IFIFECk»AY(I).GT.3C.O)25,ie 

TCHWY(K)=TCFWY(K)+HEOWAY(I) 

TCN(K)=TON(K)+ICN( I ) 

TCFF (K  )=TCFF (K ) ♦ ICFF(  I ) 

GO  TO  25 
M=1 

IF!  TCFWY!  K)  . EG.O)  TCFWiY  ( K)  = 1,0 
RACMK)=TCN  (K)/TCFWY(K) 

RACFF  !K)  =TGFF ! K ) /TCFWY! K ) 

WR  ITE!22, 24) INT !K)  ,RACN( K) ,RACFF( K) 
FCRMAT(I4,1>,F6.2,1X,F6.2) 

IF! M.EC.l )GC  TO  30 

K=K+  1 

GO  TC  15 

CCNT  IME 

END 
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APPENDIX  11 
PROGRAM  COMB  IN 


Program  COMBIN  combines  all  station  data  into  an  input  matrix. 


D I^E^SI  CN  P/CN  (20  ,28  ) »KACFF(  20,26  ),  INT  (5C  ) 

DCUBLE  PRECISICN  FILEIN 

OPEN  (UNIT=2  2,CEV  ICE=*CSK»,  ACCESS=  'SEOaUT*  , F I LE=  • PA  S CU  T.  ' ) 
OPEN  (LMT=22  ,CEV  ICE=*  CSK  S ACCESS=  'SEQOUT  • , F I L E=  ' PA  S 5 I N.  '1 
DATA(  INK  K)  ,K  = 1 ,26)  /700  , 730 , 6 00 , 83  C , 9 00  ,93  0 , 1 000 , 1 030  , 

2 1100,1130,1200,1230,1300, 1330, 1^00, 1A3C,15CC, 1530, 

2 16 CC ,1630 ,1700,1730,1800, 1830  ,1900/ 

WRIT£(5,  1 ) 

1 FCRf'AK*  NIPBER  CF  STAT1CNS  = ',  $) 

ACCEPT  2,  MS 

2 FCRMA7(I2) 

DC  50  1=1, MS 
TYPE  6 

6 FCRMAK  MNPL7  FILE  = • ,S  ) 

ACCEPT  7,  FILEIN 

7 FORMAT(AIO) 

OPEN  (LNIT=2  1, CEV  ICE= • DS K • , ACCE S S=  'SEQIN',FILE  = FILEIN) 

DC  25  J=l,26 

REAC(21,1C,EN0=12 )PACN(I,J) ,RACFF( I,J ) 

10  FCRMA7(5X,F6.2,1>,F6.2) 

GC  TC  25 

12  RAON(I,J)=0 

RAOFF  ( I,  J )=0 
25  CCNTIKE 

CL0SE(UNIT=21) 

50  CONTINUE 

DC  6C  J=l,26 

WRIT  E (22,52  ) INT  {J  ),  ( RACN  ( I , J ) , I = 1 ,MS ) 

52  FCRMAT(I4,1X,26(F6.2) ) 

60  CONTINUE 

DO  70  J = 1 ,26 

WRITE  (23,  52)  1 NT  (J)  , (RAOFF  { I,  J ),  1=1,  MS  ) 

70  CONTINUE 

ENC 
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APPENDIX  12 
PROGRAM  ESTREG 


Program  ESTREG  estimates  the  missing  values  by  regression. 


DOUBLE  PRECISION  FI LEIN , FILEOO 

DIHENSION  XBAE  (35)  , SI D ( 35)  , D (3 5)  , RI  (3  5)  , I Sf>  7E  (3  5)  , B ( 35)  , 

2 SB  (3  5)  ,T  (35)  (35)  ,RX  (1500)  ,E  (1000)  ,ANS(10)  , RATE  (2 5,1  1)  , 

2 PRATE  (275)  ,?.ES  (2  75)  ,0 PATE (27 5)  ,EPATE{275)  ,Y  (3  5,275)  , 

2 Y2  (35,27  5)  ,TVP.R(25,  2 75)  ,SVAE(11,275)  ,P(35)  ,X1  (10000) 

2 ,INT  (2  5)  ,TVAE1  (2  5,2  75)  ,SVAE1  (11,275)  ,B1  (35)  ,MK(35) 

2 , AL(275),  BE  (275)  , STRES  (275) 

1000 

WRITE (5, 1000) 

FORMAT  (‘INPOT  FILE  = « , $) 
ACCEPT  2000,  FIISIN 

2000 

FORMAT  (A10) 
WRITE(5,3000) 

3000 

FORMAT  (' OUTPUT  FI1E=',$) 
ACCEPT  2000,  FILEOO 
WRITE (5,4000) 

4000 

FORMAT (* NUMBER  OF  STATIONS=  ',$) 
ACCEPT  5000,  MS 

5000 

FORMAT  (12) 

6000 

WRITE (5,6000) 

FORMAT  (•  NUHBEF  OF  TIME  PEPIODS=  ',$) 

ACCEPT  5000,  MT 

MK=MT4-MS-1 

ML=MK-1 

MC=HT*MS 

OPEN (UNIT=21 ,DEVICE=*  DSK» , ACCESS= • SEQIN’ , FILE=FIIEIN) 
OPEN (0NIT=22 ,DEVICE=*  DS  K* , ACCESS=  * SEQOOT ' ,FILE=  FILEOO) 

4 

5 

DO  5 1=1. MT 

READ  (21 ,4)  I NT  (I)  , (RATE(I,J)  ,J=1,MS) 
FORMAT  (14,  IX,  1 1 (F6.  2)  ) 

CONTINUE 

r 

ESTABLISH  INDEPENDENT  VARIABLES  Y 
M=0 

10 

12 

DO  12  1=1, MI 
DO  10  J=1,MS 

IF  (RATE  (I,  J)  . EO.  0)  30  TO  10 
M=H  + 1 

TVAF.  (I,M)  =1 
37A?  (J,M)  =1 

PRATE (M) = A LOG (RATE  (I, J)  ) 

CONTINUE 

CONTINUE 

13 

DO  20  1=1, M 
y (1  ,L)  =PRATE(I) 

DO  13  I=1,MT-1 
Y (1  + 1, L)  =T7AF  (T,L) 

CONTINUE 

IF  (TV?  P (MT,  L)  . NE.  1)  GO  TO  15 
DO  14  I=1,MT-1 
y (1  + 1 , L)  =-i 

14 

15 

CONTINUE 
DO  16  J=1,MS-1 
y (MT  + J.L)  =SVAH  (J,U 

16 

CONTINUE 
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PROGRAM  ESTREG  (CONT.) 


17 

20 

C 

134 

135 

C 

C 

165 

250 

252 

260 

265 

270 

275 

280 

290 

295 

300 

310 


IF  (SVAE  (!1S,L)  . NE.  1)  GO  TO  20 

DO  17  J=1,HS-1 

Y (nT+J,L) =-1 

CONTINUE 

CONTINUE 

CONVERT^!  ro  DESIRED  FORHAT  OF  REG.  SUBROUTINE 
K=0 

DO  135  1=1, «K 
DO  134  J=1,M 
K = K + 1 

XI  (K)  =y  (I,  J) 

continue 

CONTINUE 

«• 

OBTAIN  CORELATION  MATPIX 

CALL  COPRE(M,MK, 1,X1,XEAE,STD.FX,R ,B1 ,D,T) 

IDENTIFY  INDEPENDENT  VARIABLE  VECTOR  ISAVR 
DO  165  1=1, ML 
ISA7E  (L)  =L*1 
CONTINUE 

CALL  ORDER (MK.F, 1 , ML, IS AVE,R X, EY) 

CALL  MIKV  (FX,ML,DET,WK,T) 

CALL  MULTR(M,MI, XBAR,STD,B1,PX, RY, ISAVE,B,SB,T, ANS) 

BRITE  (22,250) 

FORMAT (18X, 39HANALYSIS  OF  VARIANCE  FOR  THE  OEGP ESSION//) 
HRITE (22,252) 

FORMAT(2X,19HSOURCS  OF  VARIATION, 9X, 2HDF, 2X, 1 3HSUH  OF  SQUARE 
2, 3X, 1 1 HMSA  N 5QUARE,6X,7HF-VALUE) 

WRITE  (22,260)  ANS  (5)  ,ANS  (4)  ,ANS  (6)  , A NS  (10) 

FORMAT (2X,26HATTHI3UTABLE  TO  R EGRESS ION, 2X , F3 . 0 , 3 X , 

2F10. 4,5X,F10. 4,3X,F1 0. 4) 

WRITE  (22,26  5)  ANS  (8)  ,ASS  (7)  ,s  NS  (9) 

FORMAT (2X,25HDEVIATION  FROM  EE GRISSION , 3 X , F3 . 0 , 3 X , F 1 0 . 4 , 5 X , 
2F10. 4,3X,F10. 4) 

WRITE  (22,270)  E NS  (1) 

FORMAT (//1X,10HTNTERCEPT=F9.3) 

ANS  (2)  =ANS  (2)  **2 
WRITE  (22, 275)  A NS  (2) 

FORMAT  (IX, 29HMUITIPLE  CORE.  COEFF.  SQUAF.  E=F6 . 5) 

WRITE  (22,230)  ANS  (3) 

format  (1  X,  27HSTANDAPD  ERROR  OF  ESII  MATE=F9. 3//) 

WRITE (22,290) 

FOFKAT (2X,3HVAPIABLE, 10  X, 1 2UREGF.  COEFF. , 1 OX, 10HSTD. 
2EFRCP,10X,8HCOMPUTED) 

WRITE (22,295) 

FORMAT  (43X, 9HOF  COEFF.  , 1 1 X ,7 HT-V A LUE//) 

DO  310  J=1,ML 

WRITE  (22,30  0)  ISA  VE  (J)  ,B  (J)  ,SB(J)  ,T  (J) 
FORMAT(6X,I2,14X,F10.4,11X,F9.4,10X,F9.4) 

CONTINUE 
DO  313  L=1,M 
ALSC=0 
BESO=0 
ALS01=0 
BES01 =0 

DO  311  K6=2.KI 

AL  (1)  =AL  (L)  +Y  (K6,L)  *B  (K8-1) 
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PROGRAM  ESTREG  (CONT.) 


RLS0=ALSQ4-5  (K6-1)  **2 
AL501  = ALS01  *3  (K6-1  ) 

CONIINUE 

DO  312  K7  = .iT-*-1  , MK 

BE  (L)  =3S(L)  +y  (K7,L)  *B  (K7-1) 

BESO=BESQ*P (KT- 1) **2 
BES01=BES01 +3 (K7-1) 

312  CONTINUE 

DS  = DS*Y  (1  ,L)  *AL  (L)  *BE  (I) 

313  CONTINUE 
ALS0=ALC9+» LE01**2 
BES9=BESQ+EES01**2 
DSTRT=  (DS**2)  /(ALSQ*PFSQ) 

SSPE  = ANS  {7J  -DSTP.T 
FSTAT=DSTAT* (ANS (8) -1> /SSPE 
WRITE (22,314) FSTAT 

314  FORMAT  (//*  THE  INTERACTION  TEST  STAI.  , F5TP.T  : • , 1 X,  FI  0.  4 ) 

C OBTAIN  ESTIMA'^ES  FOP.  MI3SIN0  VALUES  AND  RESIDUALS 

KL  = 0 

DO  320  1=1, MT 
DO  315  J=1,MS 
KL=KL+1 
TVAP1  (I,KL)  =1 
SVAF1  (J,KL)  =1 
ORATE(KL) =PATE  (I,J) 

315  CONTINUE 

320  CONTINUE 

DO  330  L1=1,MC 
DO  325  I=1,MT-1 
Y2  (I,L1)  =TVAR1  (I,L1 ) 

325  CONTINUE 

IF  (TVAF.1  {MT,L1)  . NE.  1)  30  TO  327 
DO  326  I=1,MT-1 
Y2 (I. LI) =-1 

326  CONTINUE 

327  DO  328  J=1,HS-1 

Y2  (aT-1  + J,L1)  = SVAE1  (J,L1) 

328  CONTINUE 

IF  (SVAE1 (3S ,L1) , NE. 1)  30  TO  330 
DO  329  J=1,MS-1 
Y2 (Mr-1+J,L1) =-1 

329  CONTINUE 

330  CONTINUE 

DO  340  L2=1,MC 
T1=0 

DO  335  K3=1,ML 
T1=TUB(K3)  *Y2  (K3,L2) 

335  CONTINUE 

ERATE  (L2)  =ANS  { 1)  +T1 

IF  (OF. ATE  (12)  .NE.O)  RES  (12)  =ALOG  (ORATE  (L2)  ) - ERATE  (L2) 

STRES  (l2)  =RES(L2)  /SQP.T  (ANS  (9)  ) 

ERATE  (L2)  =EXP  ( EP, AT E (L 2)  ) 

340  CONTINUE 

WRITE  (22,350) 

350  FORMAT (//,' OHIO. PATE' , 1 X, ' FITTED  R AT E ' , 1 X , ' EES IDU AL ' , 

2 IX,  'STAND. RES. ») 

DO  352  L2=1,nC 

WRITE  (22,351)  ORATE  (L2)  , SPATE  (L 2)  ,RES  (L2)  ,STPE3  (L2) 
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PROGRAM  ESTREG  (GONT.) 


351  FORMAT  (2x,P6.2,4X.F6. 2,4X,F6. 2,5X,F6. 2) 

352  CONTINUE 

DO  353  L2=1,HC 

IF(ORATE(L2)  .NE.O)  EP.ATE  (L2)  =OEATE  (L2) 

353  CONTINUE 
1=1 

L2  = 1 

HPITE  (22,356) 

356  FORMAT (///,• ESTIMATED  UNLOADING  RATE  - ALL  RIVERSIDE  TRAINS',//, 

2 'TIME' ,2X, 'KENMORE* , IX, ' AUDIT. ',2 X, 'COPLEY* , 2 X, * ARLIB3* , 

2 2I,*30TLST.' ,1X,'PK.ST. ' ,2 X , * GOV. CTR ' , 1 X, ' H A IBKT » , 2X, 

2 'NOR. ST A* ,1X,'SC.PK. ' , 2x , • LEG HHERE ' ) 

354  WRITE  (22,355)  INT  (I)  , (ERATE  (L2^J)  ,J=0,HS-1) 

L2=L2*HS 

1 = 1 + 1 

IF  (L2.LE. HC) GO  TO  354 

355  FORMAT (I4,2X, 1 1 (F6. 2,2X) ) 

400  END 

SUBROUTINE  DATA 

RETURN 

END 
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APPENDIX  13 
ANOVA  RESULTS 


RIVOFF.MTY;  Model  (3]  analysis  of  variance  for  the  regression  CTransformed  data) 


SOURCE  OF  variation 

DF 

SUM  OF  SQUARE 

MEAN  SQUARE 

r-VALUE 

attributable  TO  REGRESSION 

31. 

146,3501 

4,7210 

10,0632 

DEVIATION  From  REGRESSION 

65, 

30,4937 

0,4691 

intercepts  0.041 

multiple  corr.  coErr,  squares. 82757 

STANUAHD  ERROR  OF  ESTIMATES  0,685 


variable 

REGR,  COEFF, 

STD, ERROR, 
OF  COEFF, 

COMPUTED 

T-VALUE 

2 

0,0421 

0,7105 

0,0593 

3 

•0,5474 

0,4024 

•1,3602 

4 

0,4067 

0,2829 

1.4374 

5 

0,6812 

0,2829 

2,4078 

6 

0,4350 

0,3115 

1,3966 

7 

•0,2319 

0.3465 

•0,6693 

8 

0,4156 

0,3154 

1,3177 

9 

•0,1246 

0,2897 

•0,4302 

10 

-0,9837 

0.3577 

•2,7498 

1 1 

•0,3746 

0.3577 

-1,0471 

12 

•0,2753 

0,4150 

•0,6633 

13 

-0,6597 

0,4010 

-1,6452 

14 

■0,3688 

0,7105 

-0,5191 

15 

-0,3714 

0,4881 

•0,7609 

16 

0.2183 

0,3097 

0,7049 

17 

•0.1597 

0,3097 

•0,5155 

18 

•0,2602 

0.3472 

-0,7493 

19 

0,3940 

0,3472 

1,1347 

20 

0.5742 

0,3456 

1.6613 

21 

0,7678 

0.3456 

2.2213 

22 

0,7146 

0.3113 

2.2955 

23 

0,4602 

0,3456 

1.3313 

24 

0,9076 

0,3113 

2,9152 

25 

0,2987 

0.3965 

0,7533 

26 

0,8251 

0,1534 

5,3781 

27 

0,0963 

0,3139 

0,3067 

28 

1,4366 

0,2711 

5,2982 

29 

0,9529 

0.2303 

4,1369 

30 

-0,7833 

0,1937 

•4.0432 

31 

0,8289 

0,1720 

4.8184 

32 

-0,7414 

0,2025 

•3,6609 

THE  interaction  TEST  STaT,  t FSTaTI  10,0026 
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ANOVA  RESULTS  (CONT.) 


NORON.MTY;  Model  (3)  analysis  of  variance  for  the  regression  (Transformed  data) 


SOnPCe  OF  variation  DF  of  square  KEAN  SQUARE  F-VALUE 
attributable  to  regression  30,  111.1235  3,7041  5,B99o 
deviation  From  regression  66,  41,4427  0,6279 


intercept*  -0,375 

multiple  CORR,  COEFF,  square*, 72836 

STANDARD  ERROR  OF  ESTIMATE*  0,792 


variable 

REGR,  COEFF, 

STD, ERROR, 

CUMpUTED 

OF  COEFF, 

t-value 

2 

-1.7435 

0,8023 

-2,1733 

3 

•1,3899 

0,8023 

•1,7325 

4 

0,3695 

0,4577 

0,8074 

5 

0,2950 

0,3573 

0.8256 

6 

0,4035 

0,4034 

1,0003 

7 

-0,2205 

0,3640 

•0 ,6057 

8 

•0,4674 

0.3681 

-1,2699 

9 

•0,3418 

0,3152 

-1.0843 

10 

0,0007 

0,3092 

0,0023 

1 1 

-0,3389 

0,3331 

-1,0175 

12 

•0,7660 

0,3984 

-1,9228 

13 

-0,4118 

0,4562 

-0,9026 

14 

•0,4933 

0,7843 

•0,6289 

15 

•0,0676 

0.4630 

•0,1460 

lb 

•0,0327 

0,4630 

•0,0707 

17 

•0,0525 

0.4630 

•0,1134 

18 

0,4997 

0,8023 

0,6229 

19 

1,3250 

0,7955 

1,6655 

20 

0,6679 

0,4604 

1,4508 

21 

0,8779 

0,3381 

2.5966 

22 

0,9597 

0.3152 

3,0446 

23 

1,1004 

0,3152 

3,4911 

24 

-0,0027 

0.3387 

-0,0081 

25 

0,9726 

0,1869 

5,2039 

26 

0,6689 

0.2695 

2.4819 

27 

1,0098 

0,2749 

3,6734 

28 

0,4497 

0,2338 

1,9238 

29 

•1,0271 

0,2553 

•4,0236 

30 

0.6594 

0,2184 

3,0192 

31 

•0,9677 

0,1979 

-4.8885 

THE  interaction 

Test  stat,  t fstati 

2,3866 
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ANOVA  RESULTS  (CONT.) 


NOROFF.MTY;  Model  (3)  analysis  of  variance  for  the  regression  (Transformed  data) 


SOURCE  OF 

variation 

DF 

SUM  OF  SQUARE 

MEAN  SQUARE 

f-value 

attributable  to  regression 

31. 

117,3748 

3.7863 

9,8920 

deviation 

From  regression 

76. 

29,0901 

0.3828 

intercept*  0.505 

MULtIPI-E  COPR,  cocrr,  S0UARE»,80139 

standard  error  of  estimate*  0,619 


variable 

REGR,  COEFF, 

STD, ERROR, 
or  COEFF, 

COMPUTED 

T-VALUE 

2 

•0,9219 

0,6248 

•1,4756 

3 

-0,2182 

0,6248 

•0.3493 

4 

1,0731 

0,3116 

3.4433 

5 

0,4788 

0,2555 

1,8745 

6 

0,6644 

0,2818 

2.3573 

7 

-0,1813 

0,2578 

-0,7034 

8 

0,0581 

0,2611 

0,2226 

9 

-0,1435 

0,2432 

•0,5901 

10  ■ 

•0.1912 

0,2251 

•0,8492 

11 

•0,5318 

0.2843 

•1,8705 

12 

•0,2032 

0,3106 

-0,6543 

13 

0,1756 

0,3560 

0,4931 

14 

0,1775 

0,4449 

0.3990 

15 

•0,4322 

0,3134 

•1,3789 

16 

0,3830 

0,3134 

1,2219 

17 

0,1668 

0,3134 

0.5322 

18 

0,1217 

0,4408 

0,2760 

19 

0,2178 

0,6198 

0,3514 

20 

•0,1526 

0,3583 

•0,4260 

21 

0,0275 

0.2611 

0.1053 

22 

•0,0830 

0,2432 

•0,3412 

23 

0,0772 

0,2432 

0,3172 

24 

•0,1637 

0.2432 

•0,6732 

25 

•1,2278 

0,1476 

•8.3207 

26 

-0,3277 

0,2146 

•1,5274 

27 

•0,3381 

0,2045 

•1.6531 

28 

•0,1479 

0,1837 

•0,8050 

29 

•0,9634 

0,2033 

•4,7392 

30 

1.6122 

0,1704 

9,4595 

31 

1,1294 

0,1549 

7,2908 

32 

0,2395 

0,1785 

1,3416 

THE  interaction 

test  STaT,  , fstati 

16,3779 

65  copies 
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